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ABSTRACT 



Simulation of inventory systems is rapidly becoming an accepted 
management tool for quantifying proposed decision rules . The Stanford 
Research Institute's Inventory Management Simulator/Evaluator is an 
excellent example of a multi-branch, multi-item simulator,, It represents 
one of the most sophisticated and comprehensive programs available today* 
In this Guide, the simulator is completely documented in such a way that 
anyone with a general knowledge of supply systems can readily follow 
the rationale of the program. For those who have a knowledge of the 
FORTRAN language and of programming techniques, a section is devoted 
to a card-by-card description of the program. Operating instructions 
and output reports are also discussed in detail. As an added feature, 
new sub-routines have been designed to produce random demand patterns 
which approximate the Normal, Exponential, Uniform, and Poisson distri=> 
butions. 
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PREFACE 



My introduction to the Stanford Research Institute’s Inventory 
Management Simulator /Evaluator occurred during a six week field assign- 
ment to that organization in the summer of 1963. At the time I reported, 
Mr. Peter H. Butterfield, Manager of the Industrial Research Division, 
suggested that I work with Mr. Michael E. Chambreau on the redesign of 
their inventory simulator. Being interested in the inventory management 
field and fresh from a year of Operations Analysis studies I gladly 
accepted Mr, Butterfield's offer. The experience has proven very valuable 
and I am sure that it will be even more so in the years ahead. In par- 
ticular, the work of that summer helped to fill that gap which always 
exists between the theoretical knowledge of the classroom and the day-to- 
day application of that knowledge in the "real" world. It certainly gave 
me a better understanding and appreciation for both the theoretical and 
practical aspects of the operations researcher's work. Learning the art 
of applying a science is somewhat akin to putting the flesh on a skele- 
ton. You begin to recognize the value of the theoretical frame-work 
upon which a problem is solved as well as the value of the solution to 
that problem. It is with the general thought of helping others to 
bridge this gap that I have written this paper and highly recommend to 
succeeding students the study of actual research reports such as these 
listed in the bibliography. 

It is my hope that this volume will serve at least three specific 
purposes: first, to provide a publication which will give a detailed 

description of the computer program and necessary operating instructions 
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such that the simulator may be utilized by the students of the Operations 
Research, the Management, and the Management Data Processing Curriculums 
in their studies and thesis works secondly, as a basic document on which 
to build a machine simulation course for Supply Corps Officers of the 
Operations Research curriculum which could be substituted for the present 
War Gaming course now required for these officers 5 thirdly, to provide a 
program which may be utilized as a model and source of basic building 
blocks with which to design more specialized programs,, 

I would like to express my thanks to Mr. Butterfield for the courtesies 
he extended me during my visit at the Institute and for permission to use 
the simulator here at the Naval Postgraduate School. Also, I am indebted 
to Mike Chambreau for his patience and instruction in programming tech^ 
niques and to Professor Rex Shudde and Commander S. W, Blandin for their 
advice and assistance on the preparation of this thesis. 
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CHAPTER 1 



INTRODUCTION 

In documenting the Stanford Research Institute’s Inventory Manage- 
ment Simulator I have intentionally directed my efforts towards providing 
a working reference for students of the Naval Postgraduate School, An 
attempt has been made to explicitly point out the assumptions and ground 
rules of the program in order to define the system being simulated. It 
is felt that by clearly defining the system^ the simulator may be used 
by students who have no training in computer language as well as by those 
who have, I have also attempted to minimize the time required for a 
student to become thoroughly familiar with the program by presenting as 
many details as possible without overly boring the reader. The mathemati- 
cal derivations of the optimal requisition and reorder level equations 
have not been reproduced in this paper. Students who are interested in 
reviewing these are referred to reports 1 and 2 of the Bibliography. I 
highly recommend to the Supply Corps Officers in the Operations Research 
curriculum the study of SRI’s Technical Memorandum No, 3« In addition 
to providing the basis of the simulator programs, it is also an excellent 
example of Operations Research at work. Chapter III which is a card-by- 
card description of the FORTRAN Language coding and program logic may be 
omitted by the student who is not familiar with the FORTRAN. 
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CHAPTER II 



THE SUPPLY SYSTEM TO BE SIMULATED 
1. General description 

The Supply System for which the simulator is designed , consists of 
a manufacturing facility, a central or "wholesale'* warehouse and a number 
of branch or "retail" warehouses. The system is intended to be a "push" 
system with a control office located at the central warehouse which 
periodically reviews at fixed intervals the system inventory positions 
and initiates action to resupply the system by placing production orders 
or resupplies a particular location by making interim shipments from the 
central warehouse stocks. There is no provision made to relocate excess 
stock at one branch to another which has a deficiency. 

The physical distribution of stock to the branches is made from the 
central warehouse which is located at the manufacturing center. Production 
orders are coordinated with the expected needs of the branches and alloca- 
tions are made to the branches at the time of completion of the manufac- 
turing process. The production quantity is physically delivered to the 
central warehouse which immediately computes the need of each branch based 
on its current inventory position and makes direct allocations as required. 
Shipments are made to each branch on a fixed schedule depending upon the 
particular requirements of a branch. There are no provisions made for 
emergency type deliveries. The central warehouse acts as a branch or 
"retail" outlet within its region as well as being a stock point for the 
entire system. 

The simulator was originally designed in connection with studies of 
products which have a common production operation up to a certain point 
and then are either modified slightly or packaged in different sizes which 
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then become distinct items of the product in the inventory,, These items 
retain a certain degree of interchangeability with respect to management 
decisions. There is, however, no provision in the simulator to substitute 
one item of stock for another to fill customer demands. However in re- 
viewing the system as a whole it is desirable to consider not only the 
items individually but also to consider the inventory status of the pro- 
duct they comprise. Examples of this type of product are lubricants, 
paints, glass, plastics, lumber, clothing, gasses, canned foods, dairy 
products, etc. Even though one particular item (or size) may not be 
readily substituted for another, there is still a management need to know 
the total quantity and/or value of the basic product for determining pro- 
duction and inventory policies as well as production quantity. For example, 
if a clothing manufacturer-distributor wanted to produce polka-dot under- 
shorts in various sizes he would need to know the total quantity of cotton 
broadcloth to produce first. (In order to accumulate product data the 
units of an item are converted to the product-unit by a size factor and 
then sum over all items related to the product.) 

The simulator can also be used for systems where this feature does 
not occur simply by considering each item as a product and setting the 
size factor to unity. 

The SRI Inventory Management Simulator /Evaluator is time oriented 
and performs the usual functions of forecasting, ordering, distributing, 
issuing, and record keeping within each time frame of the program. In 
the present form, the simulator has fifty-two time cycles representing the 
fifty-two weeks of a year. Within each time frame all the normal inven- 
tory functions are carried out in a set sequence which was arbitrarily 
chosen. This arbitrary sequence of functional operations does not bias 
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the system since over a long run the system performance will be indepen- 
dent of the actual functional sequence used within a time frame . 

The data for one product and all of its related items are processed 
and massaged completely before beginning simulation of the next product. 
Hence, there is also a cycling cn the various products involved in the 
simulation to produce the overall system data which is desired in the 
summary reports. 



2. Optimal values for system control rules 



The following variables are designed as control variables and form 
a set of parameters which will provide definite policy rules for the 
systems 



RIJK(J,K) 



XIJK( J ,K) 



XIJ(J) 



XI 



XIJK(J,1) 



The requisition objective level for the jth item of & 
product at the kth location. The maximum level in item 
units to which the inventory position is raised at. pro- 
duction time. 

The reorder level for the jth item of a product at the 
kth location. The level in item units below which an 
interim branch order is placed. 

The system-wide reorder level for the jth item of a 
product. The level in item units below which a new pro- 
duction order is placed. 

The system wide reorder level for a product. The level 
in product units below which a new production order is 
placed. 

The central warehouse reorder level for the jth item of 
a product. The level in item units below which a new 
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QI 



production order is placed. 

The minimum economic production order quantity (EOQ) 



Since the overall objective of the supply system is to minimise 
its total cost of operation optimal settings of the control variables 
could be obtained by determining each of the cost components (i„e. 5 
holding cost, production ordering cost, and shortage cost) as a function 
of the control variables and then minimizing with respect to all cf the 
variables. However, the dual functioning of the warehouse as a system 
stock point for interim orders and as a branch outlet unduly complicates 
the mathematics. To avoid this, the warehouse is considered to be an 
independent branch with demand rate equal to the sum of the direct demand 
rate plus the demand rate due to interim branch orders. This reduces our 
two level system to a one level system with independent branches. The 
optimal values for the control variables are then computed according to 

r 1 

formulas developed in t 2. . This approximation seems reasonable since 
the interim branch orders are not expected to be too great. 

Economic Production Order Quantity 

The optimal value for the economic order quantity is given in 





The "Qw" is the classical or "Wilson" economic order quantity which 



balances the inventory holding costs against the production ordering 



cost. The "f" is a correction factor to account for the demands during 



1 



Wilson, R. H. A New Method of Stock Control. Harvard Business Review, 
Vol. 5, 1926=1927 pp. 197=205. 
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leadtime and the variance of the demand and lead-time distributions. If 
the value of "Q" so determine approaches the shelf life limitations then 
a maximum limit is taken as the demands during a shelf -life period less 
one and a half times the demands during a production leadtime 0 This is 
a rather arbitrary setting for the maximum value . No minimum value is 
set. 

Product Reorder Level 

The optimal reorder level for the product is taken as the difference 
between the optimal requisition level for the product and the economic 
order quantity. The optimal requisition level for the product is the 
sum of the optimal levels for the branches. Since reviews of the invent 
tory positions are made at discrete intervals, the time at which the 
reorder levels are reached will fall randomly within the review interval. 
Then the expected demands which occur after the reorder level is actually 
reached and before this fact is recognized will be the demand rate times 
one-half the review period. Likewise, since the average number of units 
demanded per requisition is not necessarily unity the inventory position 
when reviewed will be below the reorder level by one-half of the average 
number of units per requisition. Hence, the reorder level is to be 
raised by one-half of the demands during a review period plus one-half 
of the expected number of units demanded per requisition to compensate 
for the discrete nature of the review time and demand quantity. 

Branch Requisition Objectives (Other Than Central Warehouse) 

The optimal requisition objective for an item-location is derived 
as a function of the demand rate, the total product demand rate, the 
variances of the demand and leadtime distributions, the production 
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leadtime 9 the holding cost and the cost of shortages, The mathematical 
details are given in 2,< A maximum and minimum level are set to guard 
against shelf life limitations and unusual values of the input para- 
meters. The maximum requisition objective level is set at the expected 
demands during a shelf life period and the minimum is set at the reorder 
level plus the demands during a production leadtime and a production 
cycle . 

The formula for the optimal requisition objective level iss 



R ijk ■ lnK 


jk (b ijk Q “ ^ ^ + C sj/^ 5c jk C hjk L jk^ //t | 


\ 

where f ijk 


/ 

is demands during leadtime corrected for the variances 
of the demand and leadtime distributions 


K jk 


is a variance factor 




is a probability factor 


C sjk 


is a cost of shortages 


k jk 


is another variance factor 


C hjk 


is the holding cost 


L jk 


is the production leadtime 


t 


is the production order cycle 



Central Warehouse Requisition Objective 

Like the requisition objective levels for the branches <, the optimal 
requisition objective for the central warehouse is determined as a func- 
tion of the demand rate, the variances of the demand and leadtimes 
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distributions, the branch reorder leadtime, the holding cost*, and the 
cost of shortages o However, in the case of the central warehouse the 
situation is complicated by the fact that the total demand is the sum 
of the direct demands and the branch interim orders , Since the later 
is a random quantity, it must be computed as an expected value , It is 
assumed that the "tail" of the distribution of the demands at any branch 
can be approximated by a negative exponential whose mean is taken to be 
that of the true distribution. Then the probability of an order occurring 
times the expected size of an interim shipment given that the need occurs 
is the unconditional expected demands on the central warehouse from a 
branch. The details of the computations of the probability of an order 
occurring and the conditional expected size of an interim order are 
given in J., , 

The cost of shortages must also be computed as an expected value 
since the cost of shortages pertaining to non-availability of stock to 
meet direct demands is not the same as that for demands arising from 
branch interim orders. For direct demands the cost of shortages is the 
profit lost. The shortage cost applied to the central warehouse for 
not being able to fill branch interim orders is computed as a function 
of the rate of change of the expected shortage quantity in a branch re- 
order leadtime with respect to the reorder leadtime „ As the reorder lead- 
time increases the probability of a branch shortage increases and hence, 
the expected shortage at the branch and resultant costs of shortage will 
increase. The actual cost to be lodged against any warehouse shortage 
is then the expected value of the direct and indirect demand shortage 
costs average with respect to the frequency of occurrence or probability 
distribution of the direct and indirect demands. 
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As with the requisition objective levels for the branches * a 
maximum and minimum levels are established for the warehouse to guard 
against shelf life limitations and unusual settings of the input para- 
meters, The maximum is taken as the total warehouse demands during a 
shelf life and the minimum is the sum of the reorder level and the total 
demands during a production order cycle and production order leadtime for 
the warehouse. 

Branch Reorder Level 

Although an optimal reorder level formula is derived in[2 0 j which 
is analogous to the optimal requisition objective . ormuia s it is not 
utilized in the simulator. The branch reorder level is taken to be the 
demands during a branch interim reorder leadtime , The simulator does 
provide a means of increasing or decreasing this level to study th© 
sensitivity of the overall system to changes in the reorder point. 

Central Warehouse and Item Production Reorder Levels 

The central warehouse production reorder level is taken to be the 
total demands on the warehouse during its production leadtime. 

The system~wide item production reorder level is the demands during 
a branch production leadtime summed over all branches including the 
central warehouse. 

Provisions are made in the simulator to adjust the reorder lev©! 
settings in order to study the effects of changes In these variables upon 
the overall system performance, 

3, Functional routines of the simulator 
This section will give a general description of each of the fun©“ 
tional routines within the simulator. Figure 1, is a block diagram of 
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these functional routines 



A, Start New Product 

The time cycle begins at this section and continues through the 
END-OF-WEEK HOUSEKEEPING section. The product p item and branch data 
cards are read into the program at this point,, All the parameters of 
the system have now been set and the simulation of a product can now 
commence, 

B, Initial Forecast of Demands 

The twelve monthly demands for each item-loeation are averaged to 
provide a mean weekly demand rate. There is some correlation introduced 
here between the predicted demands and the actual demands that the branches 
will experience later since both the initial predicted demands and the 
actual demands are based on the monthly demands. However s it is doubted 
that this correlation is sufficient to appreciably Influence the system 
performance except under rather erratic demand patterns, 

C, Initialize Order Queues 

The Order Queues are initialized to ensure that any data left in the 
queues from the previous product simulation is not erroneously introduced 
into the present product simulation, 

D, Compute Optimal Base Rules 

The computation of the optimal base rules Is really the heart of 
the simulation program. The decision rules or policies developed here 
will in effect make or break the system. In the present form of the simu-= 
lator the optimal requisition objectives and reorder levels are those 
published in [l,J . The subroutine "RXCOMP” 1 computes these levels based 
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on the initial mean demand rates of the item-locations, These levels 
are computed first in item units and then converted to weeks of supply 0 
In the case of the central warehouses, its optimal requisition objective 
and reorder level is given as so many weeks of total item demand (all 
branches ) 0 Likewise , the system item reorder level is given as so 
many weeks of total item demand, The optimal requisition objective for 
all branches other than the central warehouse is given as so many weeks 
of demand at the individual branch,, The product requisition level is 
taken as the item-location requisition objectives summed over all items 
and all locations a The optimal or economic order quantity (EOQ) is ex- 
pressed as the ratio of the initial economic order quantity to the square 
root of the total product demand rate, (See figure 2, for a block dia- 
gram of this subroutine). 

By expressing the initial requisition objectives, reorder levels 

and economic order quantities in these ways, they can be corrected or 

updated with the subsequent demand experience. For example, by express^ 

ing the requisition objectives for the central warehouse in weeks of 

total item demand, the corresponding number of physical units to use 

for the "current" week's requisition objective will be the "optimal" 

coverage period times the predicted item demand rat© for the current 

week. Hence, as the predicted demand rate changes so will the requisi- 

tion objective. The same principle applies to all the requisition ob- 

jectives and reorder levels. In updating the economic order quantity 

the square root of the product demand rate is used since the demand rate 

1 

enters the Wilson's EOQ formula as the square root, 

Wilson, R, H., A New Method of Stock Control, Harvard Business Review, 
Vol, 5, 1926-1927, pp, 197-205. 
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E, Predict Demand and Generate Actual Demands 

The subroutine PREDICT is used to forecast the weekly demand rates 
for each item-location, each item, and each product,, The forecasting 
scheme is a nine week moving average. There is no trend correction (or 
smoothing) applied to the forecast. The item demand rate is the sum of 
the direct demands (does not Include branch interim orders on the central 
warehouse) over all locations and the product rate is the sum of all item 
rates converted to product units. 

Because the input data is expressed in months a complicated indexing 
system is necessary to select the most recent "nine week period" on which 
to make the forecast, (See the detailed description of this selection 
process given in the card write up section for card 137) , 

The subroutine GENDEM scales the monthly demand data down to a 
weekly rate. The indexing used here to select the appropriate month 
for scaling approximates the calendar month-to-veek relationship as close 
as possible, (See the detailed description of this subroutine given in 
the card write up section for cards 913“939)<> 

Figure 3. is a block diagram showing the relationship of these sub- 
routines to the main program, 

F, Compute R and X Levels 

At the beginning of each time period in the simulation of a products, 
the requisition objectives and reorder levels are updated in accordance 
with the demands predicted for the period. Hence , these policy levels , 
as expressed in physical units, are changing each week but the "optimal" 
period of coverage remains constant. Under the SRI rules these "optimal" 
coverage factors are the outputs of the RXCOMP subroutine while under the 
alternate rules the coverage periods are given as the branch and plant 
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(or central warehouse) goals . It should be noted that under the later 
rules the reorder levels are taken as a percentage of the requisition 
objective. This percentage can be arbitrarily set by changing the valu* 
of the variable RGX on card 31. 

G. SRI Production Ordering Procedure 

The decision whether or not to place an order for new production 
is made on the basis of the system-wide inventory position of the product 
and/or of the item, and the inventory position of the item at the Central 
Warehouse o Each product and each item within a product classification 
has an assigned "optimal" system-wide reorder level. In addition, the 
central warehouse has an assigned reorder level for each Item In its 
local inventory. Each branch also has an assigned reorder level but the 
fact that a branch (other than the central warehouse) has fallen below 
its item reorder level will not trigger an order for new production. It 
will, however, trigger an order for an interim branch requisition on the 
central warehouse stock. When any one of the three inventory positions 
(i.e., the system-wide product level, the system-wide item level, or the 
central warehouse item level) falls below its reorder level a decision 
is made to place an order for new production. The system-wide inventory 
position of the item is reviewed to prevent Imbalances between Items of 
the same product since as mentioned before no provision is made to sub- 
stitute one item for another. The item inventory position at the central 
warehouse is checked to insure that any imbalance between branches (with 
respect to the same item) does not cause eventual shortages at the branches 
with the deficiency. Since there is no redistribution between branches, 
the central warehouse must maintain its stock levels even though there 
might be sufficient stock of the product and/or the item in the system. 
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Once the decision is made to place a production order it becomes 
necessary to determine what quantity to produce. Since the underlying 
philosophy of the supply system is to anticipate branch needs as accurately 
as possible and to automatically resupply them with as few interim ship- 
ments as possible the production allocations are based on the latest 
possible inventory position. Hence a deferred allocation procedure is 
utilized. The initial production order quantity is taken as the difference 
between the branch requisition objective and inventory position level* 
summed over all branches and all related items of a product ir* product- 
units, The quantity is then adjusted to the next higher production batch 
size (you might manufacture 10,100 yards of broadcloth,, but not 10,001), 

H, Common Production Ordering Section 

After determining the production order quantity of the product an 
initial (or tentative) allocation is made to the related items and an 
item-quantity allocation is also made to the central warehouse. These 
allocations are made to record the "dues from production" on the central 
warehouse records and the system-wide product and item records. This is 
necessary in order to prevent production orders from being placed again 
at the next review based on the same needs that trigger the current pro- 
duction order. It should be noted that item-quantity allocations and 
central warehouse allocations are made to the nearest packaged lot (e s g OJ 
case lots). The production order quantity, delivery date, allocations 
to item-quantities, and allocations to the central warehouse are then 
stored in a queue to simulate the production-run leadtime delay. Pro- 
duction order "dues" are not established at the branches but only for the 
central warehouse and the system-wide item and product records. This is 
equivalent to assuming that production order allocations are received at 
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the branches prior to making the allocations for the next production 
run. 

I. Fill Production Orders 

When the delivery date for a new production order comes due 9 the 
initial allocations to item~quantities and to the central warehouse ar© 
canceled and the system item needs are recomputed based on the current 
inventory positions of the branches. At this point in time the production^ 
run component of the production order leadtime has already expired and only 
that fraction due to the waiting and shipping components remain. Since 
the requisition objective was determined as a function of the total lead= 
time* the branches would be over supplied if the deferred allocations 
were based on the optimal requisition quantity. To avoid this situation,, 
the requisition objective is adjusted downward by the amount of the average 
demands during a production~run leadtime. In determining the needs at the 
central warehouse both its direct demands and branch reorder demands are 
considered. Since the quantity of the basic product produced may not be 
the same as the sum of the current item needs the item allocations are 
taken as the current needs times a proportionality factor which is the 
ratio of the quantity produced to the sum of the current item needs. In 
the simulator this quantity is then rounded to the nearest packaged lot 
(case). These quantities are then delivered to the central warehouse 
whose records are updated. Immediately afterwards 9 allocations are mad© 
to each item^location based on the adjusted requisition objectives and 
current inventory positions and placed into the branch order queue to 
await delivery at the appropriate time (see figure 5). The allocation 
to the warehouse is in effect the residual quantity not required to fill 
the branch needs. It is tacitly assumed in the simulator that th© 
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quantity in stock in the warehouse plus the quantity needs of the 
warehouse are sufficient to cover any increase in total branch needs due 
to scaling and rounding. This assumption seems valid except In the in- 
stance the warehouse does not stock the i*cem and lias had no demand and 
even in this case the probability of exceeding the production quantity 
to any significant extend would seem very small. 

J. Alternate Production Ordering Procedure 

Under this procedure production orders are placed at each scheduled 
review period. Ths production order quantity is based solely on the 
needs of the central warehouse since under this concept all branch needs 
are filled by placing orders on the central warehouse. No direct allo- 
cations are made to the branches. Hence;, this system can be thought of 
as a "pull" system. A minimum production order quantity is set to ensure 
economic plant operations. 

K. Fill Branch Orders 

This section simulates the delivery of both branch interim orders 
and production allocations at the branches. The branch order queue is 
screened to see if the delivery date of any order is due and,, if so, the 
on-hand, on-order, and in-transit records are updated accordingly. 

Figure 4« is a block diagram of this section of the program. 

L. Fill Demands on Branches 

In this section the actual weekly demands are presented to the 
branches for issue. If the on-hand balance is not sufficient to cover 
the weekly demand, stock may be issued in partial fulfillment of the re- 
quest or the demands may be canceled. If the stock out policy selector 
(IC0N3) is set at 1, the requests are partially filled, if it is set at 
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2j> the request is canceled 0 If the request is only partially satisfied 
the balance is canceled since backorders are not considered in the 
present form of the simulator „ The tallies representing the item-location 
on-hand balance s inventory position,, satisfied demands s and total demands 
received to date are updated as well as the system^, item and product in- 
ventory positions. 

M. SRI Interim Branch Ordering 

Although the system is designed to resupply the branches by alloca- 
tions from the production orders in sufficient time to meet their expected 
needs j, unusual demands may reduce a branch's item level sufficiently to 
warrant an interim resupply order. Consequently at each review time the 
item-location inventory positions are screened against their reorder 
points. If the inventory position is at or below this levels a decision 
is made to resupply the branch with the particular item(s) from the 
central warehouse. 

After making the decision to place an interim branch order we must 
determine the quantity to requisition. As in determining the branch need 
for the deferred allocation we must adjust the optimal requisition level 
downward by the average demands during the period between the placing of 
the most recently delivered production order and the present time. In so 
doing we are in effect compensating for only those demands which are a- 
bove the normal rate and hence s are allowing the system procedures to 
function with the least interference. The difference between this adjusted 
requisition level and the current inventory position is the •'optimal 8 ' in- 
terim order quantity. However s in some instances^ depending on the rela- 
tive lengths of a production order cycle and the production lead time 9 
this quantity may be too low and a minimum quantity must be set. 



17 



In determining the minimum value of the interim order quantity when 
the production leadtime is short compared to the production order cycle s 
it may be assumed that the need for reordering will occur within a pro- 
duction leadtime , (i,e,) after an order for new production has been 
placed but prior to its receipt at the brancho It may also be assumed 
that this need will occur within the latter half of the production lead- 
time, Hence, if the branch inventory position is raised to the reorder 
level plus one half of the expected demands during a production leadtime 
this quantity should suffice until the next production allocation is 
received. 

On the other hand, if the production leadtime is long compared to a 
production cycle (period between placing orders for same product) the 
deliveries will be more frequent with perhaps smaller quantities and the 
need to order will occur most likely within a production cycle and pre- 
sumably within the last half of that cycle if the system is performing 
fairly close to expectations. Hence, in this case the minimum interim 
order quantity should be the difference between the reorder level plus 
the expected demands in one half a production cycle less the actual in- 
ventory position. 

The actual interim branch order is then the maximum of the "optimal" 
quantity and the smaller of the minimum order quantities. In the simula- 
tor this order quantity is then taken as the next highest case or pallet 
lot depending upon whether the additional holding costs are greater than 
the additional handling cost incident to processing broken pallet loads. 

Once the order quantity is finalized, the simulator updates the 
warehouse records and stores the order with its delivery date in the 
Branch Order Queue (see figure 5,)» At the appropriate time the order 
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quantity is removed from the queue and the branch records are updated 
accordingly* 

N. Alternate Branch Ordering Procedure 

Since branch orders are only placed at the scheduled review times 
the simulator first checks to see if a review will be made in this time 
frame. When a review is to be conducted „ the stock status of the central 
warehouse is first cheeked to ensure availability of stock. If the cen- 
tral warehouse inventory position is below its reorder point the require** 
ment at each branch is cut in half to supply some stock to all the branches 
which require replenishment. It is implicitly assumed in both instances 
that the total quantity supplied to the branches does not exceed tne 
warehouse's on-hand balance. It should be noted that there is no provi- 
sion made for placing branch orders at times other than at the scheduled 
review date. For example, if a zero balance is reached no action is taken 
until the next review regardless of how long that interval might be. If 
the review interval is one time frame unit there are no difficulties. 
However, if the review period was several time frames in length and if 
performance of the system under the alternate rules were to be compared 
to its performance under the SRI or any other set of rules this limitation 
might be significant, 

O. End -of -Week Housekeeping 

A miscellaneous section where the average inventory on hand, the 
average value of the inventory in-transit and the number of weeks of 
stockout tallies are updated. There is also an optional print statement 
to produce an inventory status report for each item-location for each 
time period. This section also marks the end of a time frame (week). 
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The program either returns to the predict section to begin the next 
time frame or continues on to print the output reports if the present 
product simulation is completed 0 
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SCHEMATIC DIAGRAM OF SIMULATOR 



PROGRAM INV5IM 

. V 

1 S TART NEW PRODUCT 
Read input cards 



V 

INITIAL FORECAST 

Compute mean demand rate over last 
twelve months 



[ INITIALIZE ORDER QUEUES 

Zeroes out all the locations 
reserved for the production 
orders and branch orders in the 
SUBROUTINES QUEUE P and QUEUE B 



v 

COMPUTE OPTIMAL DECISION RULES 

Call subroutine RXCOMP to compute 
opt imal coverage periods 



__Y. 



GENERATE ACTUAL DEMANDS 

Scale monthly demands to weekly rate 



i 

I 

t 



. I 

FREDICT DEMANDS 



Forecast demands for the period 
! using a nine w eek moving av erage 

t V 

COMPUTE R AND X LEVELS 



Update optimal rules with current 
! demand rates 

" “1 



Figure 1 



\j 

i 

FILL PRODUCTION ORDERS 

Simulates delivery of the item from 
the plant to the warehouse and makes 
final allocations to items and to 
branches. Places allocations in the 
branch order queue. 



V 

FILL BRANCH ORDERS 



Simulates the receiving of both 
interim orders and production 
allocations at the branches 



FILL DEMAND ON BRANCHES 

Simulates issues or sales of items 
at the individual branches and up- 
dates inventory records 

V 

SRI INTERIM ERANCH ORDERING 

i 

Determines whether a branch needs 
an interim replenishment and if so 
the appropriate quantity to order 
from the warehouse 



y 

SRI PRODUCTION ORDERING 



Decides whether or not to order and 
computes order quantity 

V 

COMMON PRODUCTION ORDERIN G 

Makes initial allocations of the 
order quantity to the central ware- 
house and to items for record pur- 
poses. Places order in the produc- 
tion orde r q ueue 

v 

END OF WEEK HOUSEKEEPING 

Miscellaneous section where certain 
inventory tallies are updated 

f ' 

Figure 1 
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V 

PRINT ITEM SUMMARY^ 

Summary of system performance by 
item-location >, item;, product with 
system tota ls 



i 

[ After last prod'uci 



V 



PRINT BRANCH SUMMAR Y 

Sutnaiary of system performance by 
branches with system totals 



PRINT BASE RULES SUMMARY 


Comparison of sysi 
data between the i 
the Alternate pol 


tern performance 
SRI rules and 
Ley rules 


\ 


/ 


STOP PROGRAM 



-> 



START NEW PRODUCT 



Figure 1 
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COMPUTE OPTIMAL BASE RULES 




Figure 2 
24 



PREDICT DEMANDS AND GENERATE ACTUAL DEMANDS 




Continues 



Figure 3, 
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RECEIPT OF ORDERS AT BRANCHES 



MAIN PROGRAM 



FILL BRANCH ORDERS Section 



Sc 

qu 



V 



^ Call Queue B(2) 



reen the 
leue again 



Orders due 



SUBR OUT I NEQUEUE (N) 

Screens the order queue to determine 
if any orders are due in this time 

frame 



V 



Updates branch inventory 
records 



V 



GO TO Statement 



E 



No orders due 



V 



FILL DEMANDS ON BRANCHES 



V 

Continues 



Figure 4 
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PLACE ORDERS IK ERANCH ORDER QUEUE 




Figure 5^ 
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CHAPTER III 





DETAIL DESCRIPTION OF THE FORTRAN LANGUAGE 
CODING OF THE PROGRAM 


Card No. 


Description 


1 


Program Name 


2-4 


COMMENT cards 


5-8 


Dimension Statement. These variables are used in mala 
program in addition to those variables that are listed 
in common. There are 25 possible branches and 20 
possible items or sizes. 


9-27 


Common Storage. Dimensioning and placing of variables 
used in main program and subroutines into common 
storage locations. 


28-37 


Start of Program. Initialize various variables and 
setting End of File flag to 1. Program stops on End 
of File flag set to 2. 


38-41 


Read in the Control, Base Rules parameters, and 
Branch data cards. These decision rules and parameters 
are now set for the entire run. If different decision 
rules and/or parameters are to be used new runs will 
have to be made separately. The variable Z is used 
as a dummy for both the obsolescent cost and the 
variable VLP1 which is no longer used. 


42 


The opportunity cost is read-in in the same format as 
the physical storage cost per unit (APB and APP) and 
must be multiplied by .1 to obtain correct value. 

This is done solely for convenience of formatting. 
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43-47 Control Data.,. Decision Rules parameters and Branch, 

data are printed out to check correctness of read -is 

and also to have this information listed on the out— 
put sheets for future identification and reference „ 
48-50 Leadtime components are read-in in "days" and are 

converted here to "weeks" ^ 

51 Call SUBROUTINE RXCOMP to initialize variables vised 

for system totals „ 

52 Call SUBROUTINE RDTP to read-in the first product 

characteristics j, i 0 e. p the base number assigned to 
the product? the shelf life expressed in months j the 
branch goal expressed in daysj the plant goal (central 
whse. stock) expressed in daysj the production lead- 
time expressed in daysj the batch size increment ex- 
pressed in product units j minimum batch size expressed 
in product units „ All other product cards are read 

by card number 105 >. 

53-54 Use a do-loop to set all of the "intransit" leadtimee 

for the branches to same value 4 This is equivalent 
to assuming branches are equidistant for shipping 
purposes from the main warehouse/plant „ 

55 Re-set the warehouse shipping leadtime to zerc* Sine® 

the warehouse is located at the manufacturing plant 
there is no shipping time involved „ 

56-68 Use a do-loop to initialize branch variables 0 

69-72 Start New Products Set the time clock to zero and 

initialize the product variables « 
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73-79 



Convert shelf life* branch goal* plant goal, and 
production leadtime to "weeks", and transfer the 
values cf batch size increment (BSI), batch size 
mean (BSM), and base number (BASE) to permanent 
storage locations c It should be noted that all 
product, item, and branch characteristic data are 
first read into the AO array and then transferred 
to permanent variable name locations . When the 
program was used by SRI this data was read-in from 
a tape rather than cards and having a standard format 
for presenting information simplified the problem of 
handling massive amounts of data, 

80 Call the subroutine RDTF for the first Item of the 

present product. Here the size code, item number, 
cartons per pallet, units per cartons, weight per 
unit, production unit cost, selling price, and size 
conversion factor (# of product units per Item unit) 
are read. All subsequent item, and branch data cards 
related to the present product are read by the ’’Call 
RDTP" on card number 105. 

81-82 Set the count for the different items (sizes). As 

each item card is read the program returns to card 
82 where the item count is increased by one, 

83 Initialize the inventory position of the Item. 

84-95 Initialize the branches' (item-location) inventory 

recordsj satisfied demands, order count, stock out 
weeks, stock out counts, broken case count, broken 
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pallet count, on order quantity, average inventory 
position, total demands and the branch flag e 


96-104 


The product number (four digits) and the size code 
(two digits) are listed as a six digit number on the 
data card 0 Cards 96, 97, and 98 are used to separate 
these two numbers and store them in permanent loca- 
tion and to scale data as necessary,. 


105-124 


On the first pass through this section of the pro- 
gram we have already read the first product card and 
the first subsidary item card 0 When we call "RDTP" 
the first time at card 105 we will read a branch data 
card which will leave a value of 3 in the A0(1) loca- 
tion 0 By testing this variable, the next two cards 
will send us to card 110 where we select out the 
proper value of k by matching the branch code array 
to the branch code number read from the data card 0 
Having now determined the proper item and branch 
numbers (j and k) we transfer the input data to perma- 
nent storage arrays 0 As a bonus operation we also 
compute the mean item-location demand rate and the 
item and product inventory positions 0 


125 


This card sents us back to card 105 where we call 
RDTP again and read the next card 0 Four sequences are 
possible at this junctures 

(1) If the next card is another branch usage 
data card we will repeat the same sequence 
as outlined above 0 
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(2) If the next card read is an item record 
card, A0(l) will be set to 2,0 and this 
time card 107 sends us back up to card 82 „ 

In this way the machine recognizes that no 
more branches are carrying the previous 
item and that a new item is now being con~ 
sidered (but still under the same product). 
The item count, NJ, is increased by 1 and 
we initialize all the item and branch 
variables for the second item 0 We are then 
back to card 105 where we call RDTP in 
search of branch usage data cards for the 
new item (size). We now follow sequence 
(1) again. 

(3) If a product card is read, A0(l) will be 
set to 1.0 and card 107 will send us to 
card 129 where we commence simulating for 
the previous product and subsidary items . 
Hence, the data for the next product re= 
mains in the AO array until the simulation 
is completed for the previous product, at 
which time the program returns to card 69o 

(4) If? when card 125 sends us to card 105? a 
blank data card is read, the RDTP subroutine 
sets A0(l) to =1.0 before returning to the 
main program. Hence, when we test A0(l) 
this time we are sent to card 108? which 
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sets End of File to 2» We are then sent 
to card 129c We have completed reading in 
all the data cards and will simulate for 
the last product, print out the summary 
statements and stop the program,, 

Initialize the Order Queues, Before beginning the 
simulation we wish to initialize the Production Order- 
Queue and the Branch Order Queue, The present routine 
sets the initial values to 0,0 so in effect we ar® 
simulating a system from its beginning,. If this is 
not desirable the program could be modified so that 
representative values could be read in to simulate a 
system in progress. 

Compute R and X Levels, Here we check to see what 
set of Decision rules will be used in controlling 
the production ordering and branch shipments. If 
IC0N1 has a value of 1 we use the SRI rules, others 
wise we use the alternate set of policy rules , 

The ’’RXCOMP" subroutine is called to obtain the 
optimal policy rules for the present product. The 
requisition objectives and reorder levels for each 
item location, the system«=wide item and product 
reorder levels and the minimum economic order quantity 
are computed within this subroutine. The requisition 
and reorder levels are expressed in weeks of supply 
rather than in item or product units so that they 
may be varied as the current demands varies. 



133-135 



"Generate Actual Demand , Predict Demand v and Compute 
X and R Levels" heading , 

136 Start clock by setting time equal to 1 0 

137 Call SUBROUTINE PREDICT to make forecast of the 

weekly demand rate for each item-location, each 
item, and for the product. Demands aro computed 
by taking a weighted average of the demands in 
those months in which the previous nine weeks 
occurred. For example , the prediction for the 
first week is based on the demands occurring in 

the nine weeks (the previous two months) immediately 
preceding the simulation period. It is assumed that 
five weeks occurred in one of the months and four 
weeks in the other. Hence, the demands in the past 
two months are weighted by five-ninths and four- 
ninths. This weighted monthly average is then 
scaled by 7/30 to obtain a weekly average. In the 
Simulator, the demands of the eleventh and twelveth 
months are also used as the demands for the two months 
prior to the simulation period. 

In order to select the appropriate months in which 
the past nine weeks occurred, a complicated indexing 
system is needed. The schedule below indicates the 
results of this system for the first three time 
frames. The J1 index provides a means of assigning 
the nine week period to the appropriate months as 
the time index moves through the 52 simulation periods . 
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The DM(i) is the demands in the i-th month c 
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139 



140 



J1 Index 


Monthly Demand 


45 


DK(ll) ^ 


"N 


46 


DM(ll) 


47 


DM(ll) 




48 


DM(ll) 


1 

1 


49 


DM(12) 


T 1 


50 


DM (12) 




51 


DM(12) 




52 


DM(12) 




53 


DM(12) ) 




54 


DM(1) 


) 


55 


DM(1) 





A 



> T 



/ 



e oc 



Call SUBROUTINE GENDEM to generate the actual weekly 



demands by scaling the monthly demands inputted from 
the data cards by 7/30 and rounding to the nearest 
whole unit. Note that the predicted demands are de- 
signed as DIJK whereas actual demands are designated 
as DWIJK. 

Branch to SRI rules or the alternate set as deter- 
mined by the value of the inventory rule selector 



(IC0N1); 



Compute the economic order quantity (EOQ) for the 
current week. The variable FLQI is the ratio of the 
optimal economic order quantity to the square root of 
the mean product demand rate. By multiplying this 
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I 



I 



i 



141 



142-143 

144-145 



146 



147 



148-150 



ratio by the current predicted product demand rate 
we are varying the EOQ to reflect the present demand 
trends. The variable FLDI is the maximum EOQ due to 
shelf life limitations. 

Compute the production order cycle as the EOQ coverage 
period plus one half of the review interval. 

Initialize the product requisition and set up a 
do-loop on items. 

Compute the requisition levels for the central ware- 
house. The variable XRIJZ represents the optimal 
coverage period (of total item demands) and is con- 
verted to physical units by multiplying by the cur- 
rent predicted demands for the item. The variable 
XDIJK represents the maximum level due to shelf life 
limitations. The minimum level is given as the re- 
order point plus the demands in a production order 
cycle and leadtime. 

The reorder level for the central warehouse is taken 
as the total demands (including interim orders) during 
a production leadtime. The variable XXI JZ is the ratio 
of the mean total demand during a leadtime to the mean 
direct demand rate. 

The system-wide item reorder level is also taken as 
the demands during a leadtime. 

Initialize the item requisition variable, establish 
a do-loop on locations and screen the branches to see 
if the item is carried. 
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151 

152-153 



154 

155 

155A-156 

157 



158-159 



160 

161-167 



The item-location reorder level is taken as the 
demands in a branch interim order leadtime. 

The optimal requisition objective for the individual 
branches are computed by a formula analogous to that 
used for the central warehouse. However , the pre- 
dicted weekly demands at the branch and its produc~ 
tion order leadtime are used rather than the total 
item demand rate c (See cards 144 and 145) 

Sum the requisitions objectives for each item over all 
branches except the central warehouse. 

Continue statement for do^loop on locations 0 
Compute the total requisition level for each item 
and each product. 

The product reorder level is the product requisition 
level less the economic order quantity* 

The variable FLPDTM represents a specific time. It 
is the time at which the most recently delivered 
production order was placed. The difference between 
this time and the present time represents that portion 
of the resupply cycle which has already expired. This 
period is used to adjust the branch requisition level 
when computing the order quantity for interim ship= 
ments. (See cards 250*, 251) 

Jump statement to skip over the following section 
which is used for the alternate policy rules. 

The requisition objectives and reorder points ara 
calculated for the central warehouse and branches 
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under the alternate policy rules. The warehouse 
levels are a function of the total item demand and 





a predetermined goal. The branch levels are func- 
tions of the local demand rate for a given item 0 


168-170 


FILL PRODUCTION ORDERS Heading 


171-172 


Screen the production order queue to determine if 
any orders have been completed by the plant c If 
there is a positive order value due at this time, 
the program will continue on to simulate receiving 
the order in the central warehouse and making 
allocations to the individual branches . 


173 


Check to see which set of policy rules are in use c 
Under the SRI rules the program will compute the 
current needs of the branches and make automatic 
shipments to satisfy these needs. Under the alternate 
rules all production orders are delivered to the ware- 
house for storage „ The branches place their orders 
on the warehouse as their reorder points are reached 0 


174-179 


When production orders are initiated a temporary 
allocation is made for each item of the product based 
on the system needs at that time, and the system item 
and product records as well as the central warehouse 9 s 
records are adjusted for the production "dues". These 
initial "dues" are cancelled here prior to making th© 
actual allocations to items,, 


180 


Compute the warehouse needs for each item of the 
product. 
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181-184 


Compute the branch needs for each item and sum over 
all branch needs to give the total allocation of 
the product to each item c 


185-186 


Compute the total system current needs in product units 
and ensure that it is a positive quantity., 


187-189 


The item needs are scaled to the actual quantity 
produced (P0BS) and rounded to the nearest case lot c 


190-196 


The appropriate inventory records are updated to re- 
flect the delivery of the new items into the central 
warehouse o 


197 


Check to see which set of control rules are in use 0 


198-204 


Compute the order quantity and delivery date for 
each item-location other than the central warehouse 0 
The order quantity is rounded to a full case lot 0 
The value of inventory in transit is also updated 
here* 


205 


The PB0 (Place Branch Orders) subroutine is called to 
update the inventory records and cost tallies for 
each item-location,, The parameter value of 2 in the 
Call indicates that this is a production allocation 
and not an interim branch order „ The PB0 subroutine 
then calls the Queue B subroutine and places the 
orders in the queue to simulate the shipping delay 0 


206-207 


After completing the first production order the pro- 
gram returns to card 171 to screen the queue for any 
other orders that might be due„ 


208-210 


FILL BRANCH ORDERS Heading 
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211-215 


Screen the branch order queue to determine if any 
orders are due for receipt at the branch. If an 
order delivery date is due the on-hand , on order 
and in-transit inventory records are adjusted 
accordingly. 


216 


This GO TO statement cycles the program back to the 
CALL QUEUE B card to have another look at the queue. 


217-219 


FILL DEMANDS ON BRANCHES heading 


220-221 


Set up do-loops on items and locations. 


222 


Check to see if Branch carries this item. If not,, 
machine jumps to the next K value and checks again. 


223-224 


Check to see if there has been any demands this week 
for the item-location. If so, the total demands to 
date tally is increased by the current demands. 


225-241 


If the on hand balance is positive we compute the new 
on hand balance (TP) and the quantity of demands 
satisfied. If the new on hand Is also positive we 
jump to statement 523 to update the branch. Item and 
product inventory records. If the new on hand balance 
TP had been negative or zero the program checks to see 
which issue policy is in use. If issues are to be 
made to a zero balance (a setting of 1) the old on 
hand balance becomes the satisfied demand (TP2) and 
the new on hand balance is zero. In this case the 
stockout count tally is also increased by 1 since only 
partial issue could be made. If on the other hand the 
policy is not to make partial Issues the old on hand 
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242 

243-245 

246-248 

249 



250 



251-252 



balance is placed in TP and the satisfied demand (TP2) 
variable is set to zero and stock records are updated 
accordingly,, 

Branch to the SRI rules or alternate policy rules 0 
SRI BRANCH ORDERING PROCEDURE heading 
Set up do-loops on items and locations and test to 
see which branches carries the item. 

Review the inventory position to determine if a branch 
interim order should be placed 0 If the inventory 
position is at or below the reorder level an order 
will be placed o 

Compute the time period between the present time and 
the time at which the most recently received produc- 
tion order was placed. 

When placing a branch interim reorder the branch in- 
ventory position is brought up to its optimal requisi- 
tion level. To do so would overstock the branch and 
interfere appreciably with the automatic production 
allocation system. The policy is to supply the 
branch with enough stock so that at the time just 
prior to receipt of the next production allocation 
the branch will be at its reorder level. Hence s the 
optimal requisition level is adjusted downward by the 
expected demands during the period discussed on card 
249. A minimum reorder quantity is taken as the re- 
order level plus the demands in half a production 
leadtime (or cycle whichever is less) minus the 
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253-254 

255 

256-257 



258-259 

260-261 



262 



263 

264 



current inventory position. 

Check to see if rules call for rounding order 
quantity to the nearest carton and round if they do. 
Check to see if the rules call for rounding to pallets. 
Compute cost of issuing broken pallet load for this 
order. The cost is the number of units in the broken 
load times the cost of handling such units (cBPC), 
Compute the extra holding costs that will be incurred 
at branch if a full pallet load is issued. 

If the holding cost is greater or equal to the broken 
pallet handling cost then the order quantity is 
taken in full case lots. If not, the order quantity 
is taken as the next higher full pallet quantity. 
Compute delivery date for the order. Delivery date 
is specified as the time that the last scheduled de- 
livery was made from central warehouse plus the 
scheduled waiting period plus the in-transit time. 
Update the value of the Inventory in~transit. 

Call subroutine Place Branch Order to simulate the 
updating of the branch and central warehouse inventory 
positions, the on order quantity, the records for the 
branch, the on hand balance for the central warehouse, 
and the branch cost records. Subroutine Queue B Is 
called from within the PBO subroutine to place the 
order and delivery date in the Queue until the 
appropriate time to simulate receiving the order at 
the branch. 
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265 

266-268 

269-284 



286-287 



288-318 



319-327 



CONTINUE Statement 

SRI PRODUCTION ORDERING PROCEDURE heading 
Check Product Inventory position against product 
trigger level* If inventory is above trigger level 
then central warehouse inventory is checked against 
its trigger level. If it is above the trigger level s 
then the item inventory position is checked against 
its reorder level. If all three are above their 
trigger levels no production order is placed. How- 
ever, if anyone of the three is below its trigger 
level the production order quantity is computed by 
summing the difference between the requisition ob- 
jective and the present inventory position over all 
branches and all items s rounded to the nearest batch 
size increment. The individual item needs for th© 
system POAS(J) and the central warehouse POAW(j) are 
retained in order to make the necessary '‘dues' 8 os 
the item and central warehouse inventory records. 

This eliminates the placing of another production 
order later on to fill the same needs 0 
The production order quantity is rounded to the next 
higher batch increment. The program now jumps to 
card 319 to complete the ordering routine. The inter- 
mediate sections pertain to the alternate policy rules. 
(See Alternate Branch Ordering and Production Ordering 
Procedures at the end of this section). 

Compute the scaling factor due to rounding to nearest 
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328-329 


batch size increment (ESI) and the order delivery 
date. Production order allocations to size and to 
the warehouse are scaled to the production quantity^ 
the warehouse’s on order and inventory position, as 
well as the product and item inventory positions, 
are updated. 

A tally is kept on the number of production runs made 
during the simulation,. Finally, the order quantity 
and delivery date are placed in the queue to simulate 
the production delay 0 


330-332 


END-OF-WEEK HOUSEKEEPING heading 


333-339 


Set up do-loops on items and locations and select 
out those locations which stock the item. The average 
inventory on hand and stock out weeks records are up-- 
dated. 


340-345 


This is an optional PRINT statement used mainly for 
checking out the simulator. It provides a weekly 
status report for each item-location showing the on 
hand balance, the on order quantity and the inventory 
position. It also provides the system inventory 
position for the product and each related item. 


346 


Average Value of inventory in transit is computed. 


347 


Check the time to see if the simulation for this 
product is complete. If it is, the program will 
print out the item and product summaries. If not, 
it will return to the Generate Demand Section, set 
the time to the next increment and go through 
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348-350 

351 

352-358 



359 

360-386 



387-429 



program again. 

Print Routine for Item Summary. 

Update the page number tally. 

Set up do-loops and select out Branches carrying 
each item and print item-location data? 

DWIJK Actual weekly demands in units 

RIJK Requisition level in units 

XI JK Trigger (reorder) level in units 

XRIJK Requisition level in weeks 

XDIJK Demands in a shelf life^ item-location 
XDIJ Ratio of total warehouse demands for an 

item to total system direct demands for 
that item 

XXIJ System item reorder level in weeks 

XRIJZ Warehouse item requisition level in weeks; 

of total system direct demands for that item 
XXRIJZ Warehouse item reorder level in weeks of 
total system direct demands for that item 
XK2IK Interim branch order leadtime 

SIJK Inventory position of item-location 

Print page heading 

Initialize product and item print variables prior to 
their computation by summarizing item-location data. 

In this section the system performance data is sum- 
marized for each item with individual branch sub- 
totals. The inventory values are given in both units 
and dollar figures. The effectiveness measures such 
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430-445 



446-450 



451-463 

464 



465-467 

468-473 



474-486 



487-510 



as the stock turn ratio and the ratio of lost demand 
to total demand are based on units rather than dollars* 
The holding cost, stockout costs, and costs of handling 
broken carton and pallet lots are also computed for 
the summary report* The print statement causes the 
item-location data to be printed* 

The inventory data is summarized by product totals in 
this section. The print statement here causes the 
item totals to be printed on the summary report* 

The product effectiveness measures (stock turn ratio* 
loss demands to total demands, and weeks per stockout) 
are computed* The print statement prints the product 
totals on the summary report. 

The overall system totals (all products) are accumu- 
lated here for the Item-Product Summary Report, 

Test to see if the simulation of all products is 
complete. If not, the program returns to the START 
NEW PRODUCT section and commences simulation of the 
next product. 

PRINT BRANCH SUMMARY AND STOP heading 
The overall system effectiveness measures are computed 
and the system totals are printed for the summary re- 
port. 

Initialize the system variables prior to recomputing 
by summarizing on items and products for the branch 
totals. 

Compute the system totals by summarizing the branch 
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totals. The individual branch totals are also 
printed within this do-loop. 


511-513 


Print the system totals. 


514 


The cost of all production set ups. 


515 


The total branch variable costs for the system. 


516 


Page number tally. 


517 


Print the page heading. 


518-519 


Branch on the control rules. If the SRI rules are 
being used, the program will print a Base Rules Sum- 
mary which compares the system performance under the 
SRI rules with those of the alternate policy set. 


520 


Print branch variable cost factors and the value of 
the average inventory in transit. 


521 


Stop the program card. 


522-581 


Formats for READ and PRINT statements. 
SUBROUTINE RXCOMP (ICON) 


582 


Subroutine defined with call parameter ICON, 


583-588 


Dimension the local variables. 


589-606 


Dimension and place in Common those variables used in 
one or more subroutines of the program. 


607-610 


Function VARDF (D,C) defined. See SRI T M -3 for 
information on equation used. It computes the rel= 
variance (ratio of the variance to the square of the 
mean) of the demand distribution. 


611-612 


Computed GO TO Statement for branching into one of the 
three subdivisions of the program. The first is merely 
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an initializing routine; the second is the computation 
of the base rules themselves; and the third is the 
recap and printing of management data taken from the 
base rules. 


613-627 


Routine for initializing system management variables. 


628-631 


COMPUTE BASE RULES heading 


632-635 


Initialize variables to be computed. 


636 


Computation of expected leadtime between reaching a 
trigger level and availability of item for distribu- 
tion or the production leadtime for central warehouse . 
Sum of half of the review interval, the communication 
leadtime, and the production-run leadtime. 


637 


Set L. C. Up do-loop for locations 


638 


Compute the mean production order leadtime (reorder- 
to-receipt) for each branch. 


639 


Compute the mean branch interim order leadtime 
(reorder-to-receipt) for each branch. It does not 
include the communication leadtime or the production- 
run leadtime. 


640 


The item-location reorder point is taken as the demands 
during a branch interim order leadtime. Here the 
branch interim reorder leadtime (FLIK) is multiplied 
by the variable FK2. For management purposes, studies 
can be made of how the reorder level Influences the 
overall performance of the system by changing the 
value of this leadtime through the multiplier FK2. 


641 


T5 is an intermediate computation in the variance of 
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leadtime. It is assumed that the times when produc- 
tion is actually needed are distributed uniformly 
between times of review, and that the times of avail- 
ability are distributed uniformly between times of 
regular (scheduled) branch shipments v If the above 

assumption is true, the variance of the reaction com- 

2 , 

ponent of leadtime is FLRI 12 and the variance of 

2 

"waiting for shipment" time is FSK(K) '12 q 
T he rel-variance (ratio of the variance to the square 
of the mean) of the production order leadtime is com- 
puted as the sum of the variances of the five indepen- 
dent components (reaction time, communication time, 
production time, waiting time, and shipping time) 
divided by the square of the mean leadtime, FLIK(K), 

The rel-variance of the production-run component is 
an input parameter and is used here to compute the 
variance of the production-run leadtime 0 The shipping 
times and the communication times are taken as nearly 
constant with zero variance, 

T1 - Intermediate computational factor used to simplify 
calculations later. 

Set up do-loop on items 0 This loop is inside the do- 
loop on locations. 

Test to see if branch carries item. 

Dealer or Branch unit cost of an item, Sura the pro- 
duction unit cost CJ(J), unit freight cost, and unit 
receiving cost. 



647 



648 



649-651 



652 



653 



654 



655 



Compute the rel-variance of the demand rate as a 
function of weekly demand rate and dealer's unit 
cost. (See [l.j pages 21 and 22 for reference,) 

DlfR - Rate of arrivals of orders for the product 
(per unit time) computed as the reciprocal of the 
rel-variance of demand. (See 2. page 25) 

FIJK represents the demands during a production 
order leadtime adjusted for the variances of the 
demand and leadtime distributions. (See [l.j page 
20). T2 is an intermediate factor in computing the 
variance multiplier (FLKIJK). 

F is the sum of all the FIJK'S converted to product 
units. It is a factor in the computation of the 
economic order quantity for the product. 

Unit holding cost per unit of time is taken as the 
sum of the opportunity cost and physical holding 
cost. (Cost of holding one unit of an item for one 
year at a given branch) . 

DICHI - Factor in computing the Wilson economic order 
quantity (EOQ). In order to determine the average 
unit holding cost (over all branches and items) a 
weighted average of branch holding cost is taken. The 
weighting factor is the ratio of the item=location de=- 
mand rate to total product demand rate. The quantity 
DICHI is this weighted average times the total product 
demands rate. 

Sum the predicted weekly demands by location to give 
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656 



657 

658-660 



661-662 



663 

664 

665 



the predicted weekly demand rate for an item. 

Sum the predicted weekly demand rate by items to 
give the product demand rate in the product=equiva= 
lent units. SJ(J) is the ratio of product units to 
item unit. 

CONTINUE statement 

This FRIKT statement is used for checkout purposes 
only. The values of the print variables are checked 
against those computed manually to verify the accuracy 
of the program. These cards may be removed after the 
program is checked out. 

Computation of Economic Order Quantity heading. The 
optimal economic order quantity is computed as a 
modified "Wilson" EOQ using the following formulas 



Qw + f 



Q = < 



, 2 1/3 

1.85(Qv f) 



for 



Qw 

T 

/ 



^ 2.25 



otherwise 



where 

Qw = V 2 AD/C^j 

A is set up cost 

D is product demand rate 

C h is average holding cost per unit per year 

f is demands during a branch production 

leadtime adjusted for variances. 

Compute the square of f. 

Compute the average holding cost over all branches. 
Compute the square of the Wilson's EOQ divided by 
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666 

667-668 



669-670 



671-673 



674 



675 



f squared. The variables CM2 and FMU are multipliers 
and would normally be set to 10 and 0*0 respectively t 
Test the value of SQI to determine which formula is 
appropriate. 

If we are in the first region of the formula , we 

2 

convert the test quantity, (Qw/f) , to the appropriate 
formula by taking the square root, adding one and 
multiplying through by f. 

If the second formula is to be used then the square 
of the Qw/f ratio as computed on card 665 is multiplied 
by the constant 1.85 and raised to the one third power® 
This value is then multiplied by the factor f to ob~ 
tain the correct formulation. 

The maximum or upper limit on the EOQ is taken to be 
the demands during the shelf-life less 1.5 times the 
expected demands during the central warehouse's lead- 
time. The minimum of the optimal and maximum values 
is taken as the final EOQ for the product. This 
quantity is rounded to the next higher unit. 

The discrete adjustment period is the expected re- 
action time (FLRl/2) plus the expected time betweer. 
orders (one half of the reciprocal of the requisition 
rate) . 

The adjusted EOQ coverage time is taken as the pro- 
duction order cycle (the EOQ divided by the demand 
rate) plus the discrete adjustment period. This 
later period accounts for the fact that demands arrive 
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676 

677-678 

679 



680 



681-685 

686-693 



694 

695 



696 



in "lumps" and also for the delay between reaching 
a reorder level and taking action to reorder , 

QIFR is the number of orders corresponding to the 
adjusted EOQ quantity. 

PRINT statement for checkout purposes only. 

FLQI is the ratio of the optimal EOQ to the square 
root of the predicted demands for the product at 
the beginning of the simulation period. This ratio 
is multiplied by the square root of the current pre- 
dicted demand at the beginning of each week in the 
main program to update the EOQ value, 

FK4 and FKL are the same multiplier. The first is 
used in the RXCOMP subroutine and the later is used 
in the main program. They are both used to vary the 
leadtime demands when computing the minimum requisition 
level. (See card 720 and 153) 

Initialize the product variables. 

Set up a do-loop on items and initialize the item 
variables which are to be computed by summing over 
locations within the do-loop on branches. Note that 
the central warehouse is included in the sum 
Test to see if the branch carries the item. 

Compute the dealer's (branch) unit cost on an item 
which included the freight and receiving cost per 
unit in addition to the production unit cost. 

Compute the rel-variance of the demand rate on an 
item location as a function of the predicted demand 
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697-700 



701 



702 



703 



704 

705-706 



707 



708 



709-711 



712 



rate and the branch unit cos. 

Compute the demands during a production le&uwjsu 
adjusted for the variances due to leadtime &;id de- 
mand distributions * 

Compute the holding cost for ea:h item at each Ljc*. 

tiOHc 

FKIJK is a variance factor used in compuH.-g d,.e 
requisitioning level for each item-location,. 

The cost of a shortage at a branch is taken as the 
profit lost when cut of stock , 

The probability of ordering., 

Compute the optimal requisition objective level fex 
each item-location, (See Chapter II for delailu of 

the formula) 

Any system and/or branch constraints are applied 
here,. In utilizing these constraint factors th® 
immediate effect is to reduce the requisition ob- 
jective o This, of course will result in a lowar 
average inventory and lower holding costs 0 The 
system constraint factor may be used with or without 
any branch constraints. 

Compare the optimal requisition level to demands 
during a shelf-life and take the minimum. 

Print the item-location decision rules and input 
data for check out purposes only, 

T5 is an intermediate factor in computing the rel- 
variance of leadtime , (See card 641) 
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713 


Compute the rel-variance of the branch interim 
order leadtime. 


714-717 


Compute the demands during a branch interim order 
leadtime adjusted for the variances of the demand 
and leadtime distributions. It should be noticed 
that the factors computed here on cards 712-717 are 
analogous to those computed for the production order 
leadtine on cards 641-642 and 697-700, These vari- 
ables are used in computing the warehouse's optimal 
requisition level. 


713 


CKZJK is an intermediate factor used in computing 
warehouse requisition level. 


719 


The branch reorder level is taken as the interim 
order leadtime demands. 


720-721 


The optimal requisition level is compared to a 
minimum value which is taken as the reorder point 
plus the demands in a production order cycle and 
lead time. 


722 


Sum the item-location requisition levels over all 
branches other than the central warehouse. 


723-725 


DZIJK is the mean of a negative exponential distribu- 
tion which is used to approximate the "tail" of the 
actual demand distribution. (See 1. f pages 29-31) 


726 


FRIJK is the probability that an interim order for 
the (J,K)th item-location will occur following the 
arrival of a production allocation and before the 
arrival of the next allocation. Qr y the 
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727-728 

729 

730 

731 

732 



733 

734 

735 

736-737 



P [d in (t+FLIK) - RIJK = XIJk] . 

Compute the expected size of an interim branch 
order given that an order occurs c 

DPR(J) is the expected size (number of units) of an 
interim order for the j=th item 0 This quantity 
represents that portion of the total demands on the 
central warehouse by interim orders from the branches* 
Compute the inter arrival time of branch interim 
orders, i.e., the number of weeks per interim order. 
Compute the probability of an order for the j~th 
item. 

Compute the optimal requisition objective level ex- 
pressed in weeks. (This card would more logically 
follow after card 722). In the main program^ a new 
optimal requisition objective level is computed each 
weeks using this optimal coverage time* The new 
requisition objective level in item units is this 
coverage factor times the current predicted demands. 
The weekly predicted demands are summed over all 
branches except the central warehouse. 

The demands in a production leadtime are summed 
over all branches except the central warehouse, 

T2 is an intermediate computation to be used in 
computing the average cost of a central warehouse 
shortage. 

CSZ(J) is the average cost of a central warehouse 
shortage of the j-th item. 
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738 

7 39 “740 

741 

742 



743 

744 

745 

746 

747 

748 

749 

750 

751 



XDIJK(J,K) is the demand in a shelf=life period for 
the j-Kth item-location* 

FRINT statement for checkout purposes only* 

CONTINUE statement for the do-loop on locations. 

Compute the total warehouse demand rate as the sum of 
the direct demand rate and the branch interim order 
demand rate. Prior to this statement DPR(J) was given 
as the number of units demanded in a production order 
cycle. Here this quantity is divided by the weeks per 
cycle (STI) to give the interim order demand rate. 
Initialize the optimal requisition objective (expressed 
as weeks of supply). 

XDIJK is the demand of an item ever a shelf life 
period at the central warehouse. This quantity is 
used as the upper limit on the requisition objective. 
XK2IK(l) is the reorder time for the central ware- 
house times a multiplier. 

Compute the total direct demands per week for an item. 
Compute the total direct demands in a production lead- 
time for an item. 

Compute the total direct demands in a production lead- 
time for a product in product equivalent units. 

Compute the total direct demands in a production lead- 
time for a product considering branches only. 

Compute the total direct demands per week for a pro- 
duct considering branches only. 

The system item-reorder level is computed as system 
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demands in a leadtime. The multiplier PK5 can be 
varied to increase or decrease this quantity, 

752 Compute the total direct demands for all items of 

one product placed on the central warehouse,* 

753 XXIJ(J) is the system item reorder level expressed 

in weeks. In the main program this coverage period 
times the current predicted demand rate is used as 
the optimal reorder level for the current time period, 

754 XDIJ(J) is the optimal reorder level for the central 

warehouse expressed in weeks of total demands for the 
item in the system. It is used in the main program 

to compute the optimal reorder level in physical units 
by multiplying it by the current item demand rate, 

755 XJ1(J) is the central warehouse reorder level for the 

j-th item taken as demands in a production leadtime, 

756 XXIJZ(J) is the central warehouse's reorder level 

expressed in weeks of total item demands in the 
system. 

757 Sum the item requisition objectives to determine the 

product requisition objective. 

758 CIJ1 is the unit cost of an item at the central ware- 

house; production cost plus delivery and receiving 
costs. The later costs are given in dollars per hun- 
dred weight and must be scaled down by ,01 to give 
dollars per pound. The weight factor , WJ(J) S is given 
in pounds per unit. 

759 Transfer the unit cost figures to their permanent 
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760-761 


storage locations. 

Compute the rel-variance of the direct demands on 
the central warehouse^ and transfer this quantity to 
a permanent storage location,. 


762 


CHJJ(J) is the holding cost per unit per week,, AZR 
is the opportunity or interest cost expressed in 
dollars per dollar of inventory per year and APP is 
the physical holding cost expressed in dollars per 
product unit per year* 


763 


VDPR(j) is the rel-variance of the total demand on 
the central warehouse for each item. 


764 


DIFRFR is the requisition rate corresponding to the 
total demands on the central warehouse for each item 
summed over all items to give requisition rate for 
the producto 


765-766 


HIINT statement for checkout purposes only. 


767 


CONTINUE statement for do-loop on items „ 


768 


FK4 is a variable multiplier used to increase or de= 
crease the leadtime-demands component of the minimum 
requisition objective c It is used in the RXCOMP sub- 
routine in determining "optimal" values of the requisi 
tion objective of Item for all branches Including the 
central warehouse „ In the main program^ however , FK1 
Is used when updating the item requisition objectives 
for the branches and FK3 is used when updating the 
item requisition objectives for the central warehouse. 
This provides a means of using a different multiplier 
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on the central warehouse than on the other branches 



769 



770 



771 



772-777 



778 



779 

780-781 

782 

783 



if so desired t 

Sum the product requisition rate at the central ware- 
house to give a system requisition rate at the central 
warehouse. 

Transfer the total of branch requisition objectives 
(except central warehouse's) to another storage loca- 
tion for later summarization for system total. 

Compute the number of branch interim reorders per 
production order cycle. 

Compute the values of little f and k. Capital K is 
used in determining the optimal requisition objective. 
Compute the expected value of the cost of a shortage 
at the central warehouse. Since shortage at the 
central warehouse may mean a failure to meet direct 
demand and/or a branch interim order the average value 
must be taken. This cost factor is further complicated 
since a failure to meet an interim order may or may not 
result in an actual shortage at the branch prior to the 
receipt of the next production allocation. 

BPIJK is a computational factor representing the prob- 
ability of ordering. 

RPJK is the value of the optimal requisition objective. 
Modify the optimal requisition objective value for 
system and warehouse constraints. 

Compare the optimal value of requisition objective to 
total warehouse demands in a shelf-life period and take 
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784 



785 



786 



787-788 



789-799 



minimum quantity. 

Compare the optimal value of requisition objective 
to the lower limit and take the maximum. The lower 
limit is given as the reorder level plus a multiple 
of the total demands during a production order cycle 
plus production leadtime period. 

Add the requisition objective of central warehouse 
to total for branches to give product requisition ob- 
jective. 

Compute the ratio of central warehouse optimal item 
requisition objective to total direct demands of the 
item. This expresses the requisition objective In 
weeks of demands rather than in units. In the main 
program, the requisition objective is updated by mul- 
tiplying this factor by the current demand rate of an 
item to give the current level in units again. 

Round off the product requisition objective to the 
next higher unit and compute the product reorder level 
as difference between the product requisition objec- 
tive and the economic order quantity for the product. 
Sum the product totals to give system totals for? 

1. SD - System Direct Demands Rates All products y 

all items s all locations. 

2. SDB- System Branch Demands Rates All products,, 

all items, all branches except the central 
warehouse. 

3. SGP- System Goal for Branches? The System 
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800 



801-803 

804-814 



Branch demand rate times the branch goal 
period or the number or product units re- 
quired in a specified period (goal), 

4. SGPP- System Goal* Plant (Central Warehouse) 0 

The system direct demand rate times the 
plant goal period or the number of product 
units required in a specified period (goal). 

5. SDLB- System Demands in a production Leadtime 

for Branches only. 

6c SDL- System Demands in a production Leadtime 
all branches and warehouse . 

7, SRP- System Requisition Objective for Branches 

only, 

8, S0RD- System Order Rate 

9, SQ- System Economic Order Quantity 

10. SQP- System Economic Order Quantity for Branches 

only, 

11. SR- System Requisition Objective 

RETURN card - returns the program to the main program 
to commence simulation run on the product and related 
items for which the base rules have just been computed. 
Print Summary Section 
Management Data 

1. SG = System Goals sum of branch and plant goals, 

2. V-2 Weeks of supply in the branches. 

3. V-3 Weeks of supply in system, 

4. V-4 Difference between number of weeks of 
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supply in the branches under the SRI rules 
and under the alternate rules. 



815 

816-819 

820-825 

826 



827 

828—846 



5. V-5 Difference between the number of weeks of 

supply in the system under the SRI rules 
and under the alternate rules,, 

6. W1 - Overall System Branch Goal in weeks. 

7. W2 - Overall System Gcal in weeks . 

8. V8 - System Requisition Objective less System 

Demands in a production leadtime. This is 
the approximate maximum system on-hand 
quantity o 

9. V9 - Approximate average on~hand quantity, 

10. V10- System Requisition Objective for Branches 

less system Demands in production lead- 
time for Branches. This is the approximate 
maximum on-hand quantity in the branches, 
PRINT statement for Base Rules Summary 
FORMAT for PRINT statement of card 815, 

FORMAT statements for the PRINT statements used for 
checkout purposes. 

END (of SUBROUTINE) card. 



SUBROUTINE PLACE BRANCH ORDER (ICON) 



Definition of Subroutine 

Dimension and place in common those subscripted 
variables which are used in other subroutines of the 
program. Note that all subscripted variables used in 
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this subroutine are in common. 


847-343 


Test to ensure that the quantity to be ordered is 
positive o If quantity is zero or negative no action 
is taken to store the "order" value and program 
continues* 


849-853 


Test the branch order quantity against the warehouse 
on-hand balance to see if order can be filled* If 
not, program returns to the mainline,. However v if 
the order can be filled the warehouse on-hand balance 
and stock-position are reduced accordingly and the 
on-order tally and stock-=position of the branch is 
increased accordingly* 


854-855 


An order-count tally is kept when using the SRI rules* 
If it is desired to main this tally for other rules „ 
change the parameter ICON and read in or set the 
parameter equal to 1* 


856-859 


The following cost and count tallies are updated for 

each branch order* 

a. The cost of freight* WJ(J) is the weight per 
unit of the J-th item (lbs* per unit) $ FRK(K) is 
the freight cost per hundred weight for the K-th 
Branch and, of course,, the ,01 is a conversion 
factor from lbs* to cwts* 

b. The cost of processing (receiving)* 

c. Broken pallet count. This count is expressed in 
the number of units by which the order quantity 
exceeds the full pallet load, UPCJ(J) is the 
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number units per case and CPPJ(J) is the number 
cases per pallets. Their product is the number 
of units per pallets. 

d. Broken case count. This count is also expressed 
in units; it is the difference between the order 
quantity and the nearest full case quantity. 

e. Cost of picking and loading. This is the cost 
incident to issuing in non-pallet load or the 
cost of breaking into a full pallet load and 
issuing broken lots. CPFP is the unit cost of 
breaking the full pallet. 


860 


After updating the necessary cost and count tallies 
we are ready to place the branch order and its de- 
livery date into the Branch Order Queue. Here we are 
simulating the delay between the actual placement of 
an order and the availability of the items at the 
Branch. The subroutine parameter , N, is set to 1 to 
route us to the appropriate section of the queuing 
routine. 


861 


End of Subroutine card. 
SUBROUTINE PREDICT 


862 


Definition of Subroutine 


863-881 


Dimension and place in common those subscripted 
variables which are used in other subroutines of 
the program. No local subscripted variables are 
used in this subroutine. 
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882-383 



882-383 


The predicted weekly demands are basically a nine 
week moving average . However * since the demand input 
data is monthly demand y a complicated indexing is 
required to select the appropriate month in which the 
most recent nine weeks have occurred and the weighting 
factor given to each. Each monthly demand rate 
selected is scaled down to a weekly rats which is 
then weighted by the number of times the monthly rats 
is included in the past nine week period „ J1 and J2 
are variable indexing limits depending on time. The 
constant values 44 and 52 are chosen so that the 
first weeks prediction will be based on the 11th and 
12th month' 3 demand 0 (See card 137) 


884-887 


Initialize total product and total item demand rate 
and set up do-loops on the items and locations „ 


888 


Test to see if the branch carries the item. If the 
branch does not carry the item no prediction of 
demands is necessary. The computer will not go 
through the prediction routine for this branch thereby 
saving on computer time. 


895-397 


If the value of L as computer on card 891 is less 
than or equal to 12 we test to see if it is zero, if 
it is we add 12$ if the value lies between 1 and 12 
inclusively we use that month's demand rate scaled 
down to a weekly value. Nine such weekly demand 
rates are accumulated in the variable X. 


898-899 


The item-location demand rate is averaged over th© 
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previous nine weeks and the item demand rate is the 
sum of these over all branches. 

Continue Statement for do-loop on branch. If the 
branch does not carry the item, the program proceeds 
to the next branch without waste of time going through 
the prediction routine for that branch. 

Sum the demands for each item to give product demand 
rate. Item demand rates must be converted from units 
to product-equivalent units. 

End card for subroutine. 

SUBROUTINE RDTP (Read Tape) 

Definition of Subroutine, All the input data for 
the product, items, and branches are read into the 
program by this subroutine. The control parameters 
are read by a separate read statement at the begin- 
ning of the program. The following data is read from 
the cards: 

Product Card 

1. The product (Base) identification number 

2. Shelf life in months 

3. Production leadtime in days 

4. Batch size increment in product units 

5. Batch size mean in product units 

6. Branch Goal (Requisition Level) in days for the 
Alternate rules. 

7. Warehouse Goal (Requisition Level) in days for 



the Alternate rules. 

Item (Size) Card 

Is Size code or item identification number 
2 C ''tem conversion factor * i v e v , number of product- 
equivalents per unit of the item (e„g uy gallons 

per quart) f 

3<, Units of the item per case 

4. Cases per pallet 

5o Unit production cost in cents 

6. Unit selling price in cents 

7„ Weight per unit (pounds per unit) times ter: 

Branch Card 

1. Branch Code Number 

2 . On-hand balance for each item 

3« Inventory position for each item 
4.. Twelve monthly demand rates 

Dimension and place in common storage all subscripted 
variables which are used in other subroutines of the 
program » 

READ statement for all three types of input cards . 

The input card all have 17 fields and use a common 
format statement The data are first placed into a 
temporary storage location (A0(K))j after program 
has returned to the main program* these data are 
transferred to their permanent location „ This de- 
vice was utilized to standardize the reading of the 
input data 0 Originally all the product* iter* and 



branch data were pre-recorded on magnetic tapes 
since the volume of data was too large to read into 
the program efficiently from cards , This procedure 
could be used again by altering this READ statement 
to a READ INPUT TAPE statement, 

908 Test to see which type of card has just been read. 

A0(l) is 1 for a product card* 2 for an item card, 
and 3 for a branch card. If a product card has just 
been read the A0(l) flag indicates that it is now 
time to simulate for the previous product. 

If an item card has just been read, a positive value 
in AC(l) will send the computer back to the main 
program where the subsequent branch cards pertaining 
to the item will be read. 

If a branch card has just been read, a positive value 
in A0(l) will again send the computer back to the 
main program where the program will check for addi- 
tional items under the same product and for branch 
data under the item. 

909 After reading the last branch card of the last item 

of the last product, the program again calls the RDRT 
subroutine which attempts to read the next card. Since 
the last input card is a blank, a zero will be read 
into the A0(l) location. The test made on this card 
will cause the value to be reset to a =1. The com- 
puter then returns to the main program where the I20F 
is set to 2 signifying that the last product is ready 
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910 

911 

912 



913 

914-932 



933 



934-935 

936 

937 



for simulation and to stop the program aftsH the 
simulation is completed 0 
Return Card 

Format for READ statement 
End of subroutine card 

SUBROUTINE GENDEM (Generate Demands) 



Definition of Subroutine 

Dimension and place in common storage all subscripted 
variables which are used in other subroutines of the 
program. 

An index number is computed to determine which monthly 
demand rate will be used to compute the actual weekly 
demand rate. The schedule works out as follows. 



T L 

(Time Frame) (index for Month) 

1 * 2 , 3 , 4,5 1 

6 , 7 , 8,9 2 

10 , 11 , 12,13 3 

14 * 15 , 16 , 17,18 4 

19,20,21,22 5 

23 * 24 * 25,26 6 

27 , 28 , 29 , 30,31 7 

etc. etc 0 



Frequency of 
Selection 

4 

4 

r 

4 

4 

5 

a 4 - 

“ o 



Set up do-loops on items and locations. 

Test to see if branch carries the item. 

Compute the actual weekly demand rats by scaling the 
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938-939 


monthly demand rate down to a weekly rate* 
Continue card for do-=loops and End of subroutine 
card* 

SUBROUTINE QUEUE B (N) (Queue Branch Orders) 


940 


Definition of subroutine 


941-945 


Dimension Local Subscripted Variables „ These local 
variables are temporary storage locations for the 
order quantities, delivery dates, item and branch 
incicis. Once computed they must be stored until 
the proper time to simulate delivery at the branch,, 
One hundred and fifty storage locations have been 
reserved for these Eranch interim orders * 


946-963 


Dimension and place in common storage all subscripted 
variables which are used in other subroutines of the 


964 


program,, 

The subroutine is divided into three sections $ the 
first places the branch order in the queue, the 
second removes it from the queue, and the third see= 
tion initializes the queue* This computed GO TO 
statement directs the program to one of three sections 
of the subroutine. If the parameter value is 1, the 
subroutine places the order quantity, delivery date, 
and indices in the Queue 5 if the parameter is 2, the 
subroutine takes the orders out of Queue and simulates 
receipt of order at the branchy if the parameter value 
is 3 , the subroutine is initialized in preparation for 
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! 






the next product. 


965 


The I index counts the number of orders in the queue 
and also identifies the order. 


966 


Test to see if the queue is filled. If it is, this 
information is printed out go advise analyst but 
program continues which in effect zeros out the order. 


967-970 


The order, delivery, date, item, and location induces 
are transferred to Queue, 


971 


RETURN card 


972-973 


PRINT statement for Queue overflow condition and 
RETURN card , 


974-979 


Test to see if there are any orders in the Queue, 

If there are no orders (1=0) then the order variable 
(ORD) is zeroed and program returns to the main pro- 
gram, If there are orders in the queue, a do-Xoop 
is set up to screen through all the delivery dates 
to see if it is time to remove any orders from the 
queue and to simulate receipt at the branch. 


930-982 


Transfer the order quantity, item, and location 
indices to main program locations. 


933-986 


Transfer the last order to queue to the position 
just vacated by the order whose delivery date was 
due. This is done in order to be able to place the 
next branch order in the last queue position and to 
conserve running time and to minimize the number of 
storage locations for the queue. 


937 


Since the order data of the last queue position has 
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now been transferred to the vacated slot s the 
order delivery date (QODD) of the last position is 

erased * 


988 


The total number of orders in the queue is reduced 
to account for the order "delivered"* 


989 


RETURN card 


990-991 


This is the initialization section of the queue 0 


992 


Format statement for Print statement used to advise 
of queue overflows * 


993 


End of Subroutine Card 

SUBROUTINE QUEUE P (M) (Queue Production Orders) 


994 


Definition of Subroutine 


995-999 


Dimension local subscripted variables. 


1000-1017 


Dimension and place in common storage those sub= 
scripted variables which are used in other subroutines 
of the program* 


1018 


The subroutine is divided into three sections % the 
first places the production order data in the queue s 
the second removes it from the queue # and the third 
section initializes the queue* This computed GO TO 
statement directs the program to the appropriate 
section as indicated by the parameter value* 


1019-1023 


Use a do^loop to screen the delivery date locations 
for an empty slot to store next production order* If 
all locations are filled,, the computer will print 
"Queue overflow" statement and return to the main 
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1024-1C29 



1030-1034 



1035-1038 



1039-1040 



1041-1043 

1044-1045 



1046-1048 



program^ the order i.3, in effects cancelled , 

Transfer the production order data to the queue 0 
Allocations to sizes and the allocations to the ware- 
house for each item are stored u The machine then 
returns to the main program*, 

This is the second section of the subroutine which 
simulates the delivery of the production orders , 

Here a do-loop is used to screen the delivery dates 
in the queue to determine if any are due, If there 
are no orders in the queue the order is zeroed and 
the computer returns to the main program. 

The production order data is transferred from the 
queue locations to main program locations. 

The delivery date location of the production order 
removed from the queue is zeroed. Since there are 
only delivery dates there is no need to maintain the 
orders in consecutive locations as in the Queue B 
routine. 

Initialize the production order delivery dates 
locations. 

Format for Print Statement which advised of queue 
overflows and End-of-Subroutine card, 

MISCELLANEOUS FUNCTIONS 

Definition of RNDNF Function, This function rounds 
the input quantity to nearest lot size^ e,g„, round- 
ing an order of so many units to the nearest case lot 
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quantity. 


1049-1051 


Definition of RNDCF Function, This function rounds 
the input quantity to the next higher lot size. 


1052-1054 


Definition of END If Function, This function rounds 
the input quantity to the next lower lot size. 


1055 


End of program card 

ALTERNATE DECISION RULES 


288-290 


Branch Ordering Procedure heading 


291-292 


Determine whether it is time to conduct a review of 
the stock records. If it is time,, X will be zero. 
If it is not review time then the program jumps to 
the next section. 


293-300 


Set up a do-loop on items to screen the reorder 
levels against the inventory positions for each item 
in the warehouse. If the reorder level is less than 
or equal to the inventory position it is assumed that 
the warehouse's stock will be sufficient to meet all 
branch needs and we then test the branch reorder 
level against their inventory position for each item. 
If the branch inventory position is below or equal to 
the reorder level we compute the order quantity as 
the difference between the branch requisition and 
branch inventory position. If the warehouse inventory 
position had been less than its reorder level indicat- 

* 

ing that it was low on stock of the item then the 
branch order is reduced by one half. 
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301 


The delivery date is computed by determining the 
last scheduled delivery date then adding the time 
period until the next schedule shipment and the in- 
transit time. 


302 


The value of the inventory intransit corresponding 
to this order is added to the previous total,, 


303 


Call the P 1 ace Branch Order subroutine to simulate 
the filling of the branch order at the central ware- 
house and updating of the stock records , The PBO sub- 
routine calls the Queue B routine to place the Branch 
order in the queue until the desired delivery date 
arrives. 


304 


CONTINUE card 


305-307 


Production Ordering Procedures heading 


308 


Initialize the Production Order Batch Size variable. 


309-312 


Compute the production quantity as the sum of the 
difference between the requisition level and inven- 
tory position at the central warehouse. Production 
orders are based solely on the item needs at the ware- 
house under these rules 0 


313-318 


Test to see if production order quantity is positive. 

If it is the value is transferred to the "ord” 
storage location and then the production order quantity 
is rounded to the next higher batch size. To ensur© 
that a minimum economic order quantity is produced s 
the order quantity is taken as the maximum of the 
rounded order quantity and a predetermined minimum 
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CHAPTER IV 



OUTPUT REPORTS 

1, Control and branch card input data 

A. Selectors 

1. ICONlo Rule Selectors A value of 1 indicates that the SRI 
decision rules are to be used to compute the requisition 
levels j> reorder levels? and economic order quantities? a 
value of 2 indicates that the alternate set of decision rules 
are to be used 0 

2. IC0N2 o Not used in the present form of the simulator 0 

3. IC0N3<> Stockout Policy Selectors A value of 1 indicates 
that demands are partially filled? i 0 «o.o issues are made 
until a zero balance is reached? a value of 2 indicates that 
demands which can not be fully satisfied are cancelled,, 

4. IC0N4o Round to Cartons Selectors A value of 1 indicates 
that Branch interim orders are not rounded to case lots? a 
value of 2 indicates that the order quantity is rounded to 
the next higher full case lot„ 

5. IC0N5„ Round to Pallet Selectors A value of 1 indicates 
that Branch interim orders are not rounded to pallet lots? 

a value of 2 indicates that the order quantity will be rounded 
to the next higher full pallet lot if it is economical to do 
so„ To make this decision the holding cost of the additional 
units is compared to the additional costs incurred due to 
breaking and handling less than full pallet loads 0 

B. Costs 

1. Opportunity cost (AZR)s The cost of inventory investment in 
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cents per dollar per year, 

2. Physical Storage costs (APB and APP) ; The cost attribu- 
table to the physical handling and sheltering of the inventory 
is given in cents per product-unit per year, A different cost 
rate may be assigned to the central warehouse than that which 
is assigned to the branches, 

3. Order Preparing cost (CPPO)s The administrative cost of 
preparing and processing a production order in cents per 
order, 

4. Picking and Loading (CPL)s The additional cost of the physical 
handling of broken cases and pallets is given in cents per 
item-unit handled in a broken case or broken pallet lot. 

Note that there is no additional cost for handling full pallet 
lots since this is taken as the standard handling procedure, 

5. Cost of Receiving (CR): The costs incurred in the actual 

receiving, handling-in, and storing of items in cents per 
hundred weight, 

C. Multipliers 

1. Demand Forecasts This multiplier may be used to increase or 
decrease the predicted demands (see card 898) and thereby 
study the effects of changes in future demand patterns, 

2. Lead Times These multipliers are used to study the effects 
of increases or decreases on the branch interim reorder lead- 
times and on the production-run component of the production 
order leadtime , (See card 76) 

3. Goals s Multipliers used to vary the branch and central 
warehouse requisitioning goals, (See cards 74 and 75) 
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4. Standard Costs Multiplier used to vary the standard unit 
production cost. (See card 102) 

D. Lead time Components 

1. Reaction (LR); The time between the reaching of a reorder 
level and the recognition of that fast and the initiation of 
action to place an order. The input is in days, 

2. Communcation (LC); The time between reaching a decision to 
place a production order and the receipt of that order by the 
production unit. The input is in days. 

3. In Transit (LS); The shipping time between the central ware- 
house and a branch. In the simulator it is assumed for sim- 
plicity that the shipping time is the same to all branches. 

No provision is made for emergency shipments. The input is 
in days. 

E. Base Rules Computation Parameters 

1. AZR, APB, APP, CR; See section B above. 

2. LR, LC , LS: See section D above. 

3. Al: The cost of setting up a production run in dollars per 

set up. 

4. A3: Not used in the present form of the simulator. 

5. MU: The system inventory constraint factor. This multiplier 

is used to reduce the "optimal" requisition objective of each 
item-location by some multiple of the demands in a branch 
production leadtime. (See cards 70? and 782). It is also 
applied to reduce the economic order quantity (see card 665). 

By varying the value of MU ths system average holding cost or 
system average on hand quantity can be set to any pre-determined 
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valua, Setting MU equal to zero is equivalent to placing no 
constraints on the system, 

6. K6s Multiplier used to increase or decrease the central 
warehouse item reorder level in the RXCOMP subroutine , (See 
card 755). 

7. K5: Multiplier used to increase or decrease the system 

item-reorder level in the RXCOMP subroutine 0 (See card 751) , 

8. K4: Multiplier used to increase or decrease the item-location 

minimum requisition objective. The minimum level is taken as 
the reorder level plus the demands during the production cycle 
plus production leadtime period. The multiplier increases 
only the later component of the minimum requisition objective. 
Hence, doubling the multiplier will not necessarily double the 
minimum requisition level. This multiplier applies to the 
central warehouse's minimum level as well as to the other 
branches, (See cards 720 and 784 in RXCOMP subroutine). It 
should be noted that this multiplier is used as a dummy vari- 
able within the program and is not an input quantity. Prior 
to computing the "optimal" requisition objective level for the 
central warehouse, the value of the multiplier FK3 is put into 
FK4 (card 768), and prior to computing the level for the other 
branches the value of FKL is transferred to FK4 (card 680), 

It is not known why the variables F1Q. and FK3 were not used 
directly in computing the minimum requisition objectives, 

9. K3i This multiplier is used in the main program (card 145) 
to update the "optimal" central warehouse requisition level. 

It is used to increase or decrease the lower limit on the 
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requisition level. See K4 for further information, 

10. K2s Multiplier used to increase or decrease the item=location 
reorder level for all locations except the central warehouse 
(see card 640, 719j> and 151) , By varying this multiplier the 
lower limit on the requisition objective is also changed since 
the minimum requisition level depends on the reorder level, 

11. KLs This multiplier is used in the main program (card 153) 
to update the "optimal'* item-location (except for central 
warehouse) requisition objective. It is used to increase or 
decrease the lower limit on the requisition level. See K4 
for further information, 

12. CM2 s Cost multiplier used to increase or decrease the pro= 
duction setup cost (card 665) <> 

13. VLPls The rel=variance (ratio of variance to the square of 
the mean of the prodiction-run component of leadtime). It is 
a dimensionless quantity. 

14. VLP3? Not used in the present form of the simulator. 

F. Branch Data 

1. Codes The activity accounting numbers assigned to the 
branches. Code 1 is assigned to the central warehouse, 

2. Cost of Shipments The freight rate in cents per hundred 
weight. The freight rate at the central warehouse is gives 
the negative value of the receiving cost in order that the 
unit cost at the warehouse will be the unit production cost, 

3. Scheduled Shipping Frequency? The period between shipments 
to each branch, 

4. Branch Inventory Constraints Multipliers? If it is desired 



81 



to limit the individual branch total average inventory to 
some specified dollar value or on~hand quantity, vary this 
parameter on a trial and error basis until the proper setting 
is obtained „ Each branch can be given a different constraint 
factor if it is so desired 0 

2, The item, product, and system summary report 
The average inventory values, sales, costs, and stock effectiveness 
data are summarized here by item, product, and system 0 Complete branch 
data is listed for each item 0 The "Holding Cost" given in the report 
includes the opportunity cost factor as well as the physical holding cost 0 
The "Order Count" is the number of interim branch orders for the item* 

The "Estd Cost" of stockouts is the profit loss due to lost sales „ The 
number of broken cartons and pallets is given in item units rather than 
in the number of cartons or pallets actually opened to fill an order „ 

The other headings on the report are self-explanatory,, 

3o The branch summary report 

The inventory values, cost data, and stock effectiveness measures 
are summarized by branch for all items and products,, 

All the data given in this report under the "Inventory Values" 
heading including the stock turn ratio are expressed in terms of dollars „ 
The "Order Count" is the number of scheduled deliveries to the branches „ 
The "Preparation and Processing" cost is the cost of making the scheduled 
deliveries to the branches „ It is the number of scheduled deliveries per 
year times the cost of preparing and processing an order (CPPO)„ The 
"Picking and Loading" cost is the sum of the costs incident to handling 
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of broken cartons and broken pallet lots and the cost of handling the 
full pallet lot portion of the orders . The "Unloading fr'ocessing" cost 
is the cost of receiving and storing materials at the branch. The 
other headings on the report are self-explanatory. 

4° The base rules summary report 

1. D: The expected demands per week for the system,, 

2. RPR/D: The requisition objective in weeks for all branches except 

the central warehouse under the SRI rules . 

3. R/D: The requisition objective in weeks for all branches including 

the central warehouse under the SRI rules . 

4. AFR/D and A/D; Same as 2 and 3 above except these are for the 
alternate policy rules . 

5. (RFR-AFR)/D and (R~A)/D; The difference between the requisition 
objective requirements under the SRI and the alternate policy rules » 

6. DL(FR): The "pipeline" or leadtime quantity in product units for 

all branches except the central warehouse 0 

7. (DL) : The "pipeline" or leadtime quantity in product units including 

the central warehouse. 

8. (ON-HND-MX)? The average maximum on-hand level for all branches 
including the central warehouse. 

9. (ON-HND-MX)PR; Same as 8 above but does not include the central 
warehouse. 

10. ( ON-HND-AVE ) s The expected on-hand inventory level for all branches 
including the central warehouse. 

11. (ON-HND-AVE) Ft s Same as 9 except the central warehouse Is not 
included. 
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12. ORD: The average number of production orders per week 0 

13. Branch Picking Costs The system=wide cost of picking and loading 
for both branch interim orders and production allocations „ 

14. Production Costs The total set~up costs „ 

15. Shortage Costs The loss of profit due to nonavailability of items 
to meet demands as they occurred „ 

16. Holding Costs The physical holding costs incurred at all branches 
including the central warehouse 0 

17. Average Value in Transits The dollar value of the average number 
of units in transit between the central warehouse and the branches » 
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AZR .205 APB .240 APP .430 CR .150 
LR 4.000 LC 1.000 LS 11.000 A1 400.000 
A3 .000 MU .000 K6 1.000 K5 1.000 
K1 1.000 K2 1.000 K3 1.000 CM2 1.000 
VLP1 .000 VLP3 1.000 
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CONTROL CARD INPUTS — RUN NUMBER 
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BRANCH DATA 



CHAPTER V 



CONTROL CARD FORMATS 

Card 1. This card contains the run number, the control selectors 

settings and the system cost factors, 

a. Columns 1-5, The letters "IEENT" are placed in these columns 
to identify the run number. If desired these columns and/or 
any columns on subsequent cards used for alphabetical identic 
fication may be left blank, 

b. Columns 6-8, These are usually left blank but may be utilized 
for any non-input data as the letters "IDENT" in columns 1-5, 

c. Columns 9-14, The run number designator is entered here. It 
may be an alphanumeric designator, 

d. Columns 15=24o The letters "SELECTORS" are placed in the 
columns, right adjusted to leave column 15 blank, 

e. Columns 25=26, 27-28, 29-30, 31=32, 33=34« The values for 
IC0N1, IC0N2, IC0N3, IC0N4, and IC0N5 (i,e,, the number 1 or 
2) are placed in these columns, right adjusted, 

f. Columns 35=40, The letters "COSTS" are inserted here, also 
right adjusted, 

g. Columns 41=44, The opportunity cost (AZR) is inputed in 
mils, right adjusted, A value of 205, for example, indicates 
an opportunity cost of 20,5 cents per dollar of investment per 
year, 

h. Columns 45-58, The physical holding cost for the branches 
(APB) is inputed in cents, right adjusted, A value of 24, for 
example, indicates a holding cost of 24 cents per product unit 
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per year 0 

i. Columns 49-52 , The physical holding cost for the central 
warehouse (APP) is inputed in cents, right adjusted, A value 

43s f° r example, indicates a holding cost at the warehouse 
of 43 cents per product unit per year, 

j. Columns 53-57, The cost of preparing and processing an order 
(CPPO) is placed here in cents , right adjusted. The value of 
75 s for example, indicates a cost of 75 cents per order 0 

k. Columns 58-62, The cost of handling broken pallets (CBPC) is 
given here in cents, right adjusted. The value of 7, for 
example, indicates a cost of 7 cents per item unit handled, 

l. Columns 63-67 0 Same as k above except for broken cases (CBCC), 

m. Columns 68=72, The cost of receiving (CR) in cents per 

hundred weight, right adjusted, 

n. Columns 73-77, The cost of handling (or picking) a full 
pallet lot (CPFP) given in cents per pallet, right adjusted. 

Card 2. This card contains the management multipliers, the reaction, 

communication, and shipping components of the production leadtime and 

the cost of a production set-up, 

a. Columns 1-11, The letters "MULTIPLIERS" is inserted her© for 
identification of the card and the subsequent input data 0 

b. Columns 12-15, 16-19, 20-23, 24-27, 28-31, 32-35, The values 

for the demand forecast (FMDF), the warehouse to branch (FMLWB), 

and the plant to warehouse leadtime (FMLP) multipliers, the 
branch and plant goal multipliers (FMFG, FMBG), and the stan- 
dard cost multiplier FMCS are placed in these columns in the 
given order. The multipliers are specified to two decimal 
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places within the programs hence, a value of 100 indicates 
a multiplier of 1.00, A setting of 1„C0 for the multipliers 
is equivalent to not using them* They should never be set at 
zero 0 

c. Columns 36=45 » Leave blank 

d. Columns 46=55 . The words "LEAL TIMES" are inserted here,, 

e. Columns 56=58, 59=61 „ The value of the reaction (FLR) and 
communication (FLC) components of leadtime are given here in 
days, right adjusted 0 

f. Columns 62=64» Leave blank 

g. Columns 65=67 0 The value of the shipping or in=transit 
component of leadtime (FLS) is given here in days, right ad= 
justed 0 

h. Columns 68=72„ The letters "Al" indicating the set=up cost 
are placed here, right adjusted „ 

i. Columns 73=75 « The value of the set-up cost is inputed here 
in dollars per set up, right adjusted 0 

Card 3. This card contains the input data for the system capacity 
constraint (FMU), the reorder level multipliers (FK6, FK5), the 
minimum requisition objective multipliers (FK1, FK2, FK3), which ar® 
referred to as "floors" in [l„] , the set up cost multiplier (CM2), 
and the rel-variance for the production-run component of leadtime 
(VLP) . 

a. Columns l-3<> The letters "MU", left adjusted, are listed her® 
for identification only„ 

b. Columns 4=8 „ The value of the system capacity constraint (FMU) 
is given here to three decimal places (i<>e 0 , a value of 8300 
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Card 



Is read in as 8,300) „ 

c. Columns 9=24o The letters "TRIGGER LEVELS" are entered her® 
right adjusted , for identification only, 

d. Columns 25=27 , 28=30, The values of FK6 and FK5 are entered 
here to one decimal place, (i 0 e o , a value of 14 is read as 
1°4)« They should never be set to zero, 

e. Columns 31=38, The letters "FLOORS" are placed here, right 
adjusted, for identification only, 

f. Columns 39=41* 42=44* 45=47, The values of FKL, FK2, FK3 
are entered in these fields in their respective order. The 
multipliers are given to one decimal place (i 0 e 0 $, a value of 
060 is read as 6,0), If these multipliers are set at 1,0 it 
is equivalent to not utilizing them. The multiplier FK2 
should never be set to zero, 

g. Columns 48=59, The letters "VARIATIONS", right adjusted, are 
placed here for identification only, 

h. Columns 60=63, The value of the set=up cost multipliers (CM2) 
is entered here to two decimal places, (i,e,, a value of 404 
is read as 4«04), 

i. Columns 64=71, The value of the rel-variance of the produce 
tion-run component of the production leadtime is entered here 
to two decimal places, right adjusted, A value of 150 is 
read as 1,50, 

4. Branch Data heading card to identify subsequent input cards, 

(Notes The READ format statement is designed to skip a card at this 

point and would mis=read all the subsequent cards if this particular 

one is omitted) 
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Card 



Card 



Card 



Card 



Card 



5. Branch Code Card;, Three columns are used for each branch 
code with the code numbers right adjusted „ Any set of code numbers 
(each not exceeding three digits) may be used, There is a limita- 
tion of 25 branch code numbers. 

6. Freight Rate Card, The cost of shipping ,FRK(K) , in cents 
per hundred weight is entered for each branch in 3 ame order as the 
branch codes are entered. Three columns are allowed for each freight 
rate cost. The limitation is also 25 entries, with the numbers right 
adjusted , 

7. Shipping Schedule Card. The scheduled shipping frequency for 
each branch, FSK(K) , is listed on this card* (i,e», the number of 
weeks between shipments from the central warehouse to the individual 
branches. Three columns are allowed for each field. The values 
are given as integers and are right adjusted. 

8. The Branch Capacity Constraints, The individual branch con- 
straint multipliers (FMUK(K)) are entered here to one decimal place. 
Three columns are allowed for each entry and numbers are right 
adjusted, A value setting of zero is equivalent to placing no 
constraints on the branches, 

9. Product Data Card 

a. Column 1. The number 1 is always placed here to indicate 
that this is a product data card, 

b. Columns 2-5, The base identification number of the product 
is entered here, 

c. Columns 6-20, Leave blank 

d. Columns 21-25, Enter the shelf-life in months* right adjusted, 

e. Columns 26-50, Leave blank 
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f. 



h. 

i. 

J. 

k. 

Card 10. 
a. 



b. 



c. 



d. 



e. 



f. 



g. 



Columns 51~55o Enter the length of a production run (FLP) in 
days, right adjusted ® 

Columns 56=60,, Enter the batch size increment (BSI) in 
product units, right adjusted c 
Columns 61=65 ® Leave blank. 

Columns 66=70® Enter the minimum economic batch size (BSM) 
in product units, right adjusted ® 

Columns 71=75® The Branch Goal (BG) is given in days, right 
adjusted. 

Columns 75=30® The Plant Goal (FG) value is given in days, 
right adjusted. 

Item Data Card 

Column 1„ The number 2 is always placed here to indicate that 
this is an item data card. 

Columns 6=9® The product base identification number is entered 
here, right adjusted. 

Columns 10=11® The item identification number is entered her©, 
right adjusted® 

Columns 15=20® The item=unit to product=units conversion 
factor (or size factor) is entered here to four decimal places® 
The value of 11057 is read as 1®1057® 

Columns 51=55® The number of units per case (UPCj(j)) of the 
item is entered here, right adjusted® 

Columns 56=60® The number of cases per pallet for the item 
(CPPJ(J)) is entered here, right adjusted® 

Columns 61=65. The unit production cost (CJ(J)) of the item 
in cents is entered here, right adjusted® 
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h„ 



i. 



0 . 

Card 11. 
a. 



b. 



c. 



d. 



e. 



Columns 66=70, The sales price of the item per unit ( P J i J ) ) 
is entered here^ right adjusted. 

Columns 71=75, The weight per unit (WJ(J)) of the item in 
pounds is entered here^ right adjusted. 

Notes All other columns are left blank 
The Item-Location Data Card 

Column I, The number 3 is always entered here to indicate 
that this is an item-location data card 0 
Columns 12=X4o The branch code (BCK(K)) is entered here,, 
right adjusted. 

Columns 15~20, The initial on hand balance in item units 
(OHJK(J p K)) is entered here;, right adjusted. 

Columns 21=25 , — 9 76=80, These are 12 fields of five 
columns each in which the monthly demands are entered in 
item units. The average long run demands based on these 
numbers will be used to make the initial forecast of demands 
and to compute the initial optimal inventory policies. They 
are also used individually to compute the actual weekly 
demands lodged at the branches. 

Notes All other columns are left blank. 
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CHAPTER VI 



PROBABILISTIC DEMAND PATTERNS 

Four subroutines have been designed to utilize probabilistic demand 
patterns for the SRI simulator. These subroutines will provide 12 random 
numbers from the appropriate distribution desired for each of the monthly 
demands. The procedures for scaling the monthly demands to arrive at the 
predicted and actual weekly demands remain the same. 

Normal, exponential, uniform, and Poisson distributions are approxi- 
mated in these subroutines. The distribution parameters are defined 
within the subroutine and can be changed to any desired value by changing 
the appropriate cards. The distributions are truncated to prevent un- 
realistic values of demand from occurring. In the case of the continuous 
functions, the random numbers are rounded to the nearest integer. Even 
with these modifications a close approximation to the re3l distributions 
has been attained which is more than accurate for the purposes intended. 
The normal distribution is truncated two standard deviations above and 
below the mean. The exponential distribution is truncated at that value 
above which only a certain desired probability (error term) occurs. The 
uniform distribution, of course, is bounded by its upper and lower end 
points which may be set to any desired values. The Poisson distribution 
is terminated at three times its mean. Table 1 gives the results of 
some test runs on these subroutines. 

In order to use the random demand subroutines the program and input 
cards must be modified as follows: 

1. Input Cards, On each branch data card leave columns 21 through 
80 blank instead of inserting the monthly demand values. 
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2. Subroutine Parameter Cards, Make new cards with the desired 
values of the distribution parameters and in the case of the 
exponential subroutine set the desired value for "ERROR", 
i.e., the probability area to be truncated (a value of „01 to 
.05 is suggested), 

3. COMMON Cards, All COMMON cards of the main program must be 
listed in the subroutine and similarly all COMMON cards of 
the subroutine must be listed in the main program, 

4. Call Cards 

a. Initialize the Random Number Generator, In order to vary 
the random numbers used in different runs the random de- 
mand subroutine may be called a few times at the beginning 



of each run by the following routine; 

ISPIN = ( ANY DESIRED FIX POINT NUMBER) 65A 
CALL UNIFORM (l) 65B 
DO 99 K = 2, ISPIN 65C 
99 CALL UNIFORM (2) 65D 



b. Call for the Monthly Demands, Insert the distribution 
call card after card number 118 j for examples 

CALL NORMAL (2) 118A 

The above procedures will introduce only one random demand pattern 
for all products and items since the distribution parameters are fixed 
for each simulation run. If it is desired to vary the parameters for 
different products, items, or locations, this may be accomplished by 
listing the parameters on either the product, item or branch data card. 
Since each distribution has two parameters columns 26-30 and 31=35 of 
the data cards may be used for this purpose. The RDTP subroutine will 
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place these values in the A0(15) and A0(l6) locations. 



Products’, To vary the parameter values for each products use the 
following routine: 

1. Insert the value of the two parameters in columns 26-30 and 
31-35 in the product data card, right adjusted. These values 
will be read as integers and must be scaled if decimal places 
are desired, 

2. Place the names of the parameters in COMMON and remove the 
parameter cards from the subroutine, 

3. Prepare cards for the main program to transfer the parameter 
values from the AO array to their permanent storage locations 
For example, if using the SUBROUTINE NORMAL, prepare the 
following cards: 



FMSAN = A0(15) 


13k 


STDDEV = A0(l6) 


73B 


STD DEV = STDDEV * .01 


130 



Insert these cards after card 73 in the main program. The 
last card is an example of scaling a value to a desired decimal 
figure and may or may not be required. 

Items: To vary the parameter values for different items under the 

same product, use the following procedures 

1. Same as (l) for products except insert parameters values on 
the item data cards instead of the product data cards, 

2. Same as (2) for products 

3. Same as (3) for products except number the transfer cards 
98A, 98B, and 98C and insert these cards after card 98 in 
the main program. 
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Items^locations s To vary the parameter values for each location 
which carries the item, use the following procedures 

1. Same as (l) for product except insert the parameter values on 
the branch data card instead of the product data card, 

2. Same as (2) for products * 

3. Same as (3) for products except number the transfer cards 
116A, 116B, 116C and insert these cards after card 116 in 
the main program,, 

SUBROUTINE UNIFORM,, Uniformally distributed random numbers are 
generated on the interval ( A j, B) by the formula,, 

U(A,B) = A + (B~A)*RN 

where RN is a U(0,l) random number produced by SUBROUTINE RAND, The 
above formula was derived by equating the integral of the uniform density 
between the lower end point and a variable upper limit to the U(0,1) 
random number generated by RAND and solving for the variable upper limit 
The variable "FLWPT" is the lower end point in the floating point mode 
and the variable "UPENDPT" is the upper end pointy also in the floating 
point mode. These variables may be set at any desired values, 

SUBROUTINE EXPONTL, Exponentially distributed random numbers are 
computed using the formula , 

Y = g Log e ( l./(l. - RN) ) 

where RN is a U(0,l) random number produced by SUBROUTINE RAND, This 
formula was also derived by the integral method. The variable ’‘THETA" 1 
is the parameter of the distribution and not the mean which is, of 
course, the reciprocal of the parameter. The variable "ERROR", used 



96 



elsewhere in the subroutine, is the value of the probability area under 
the density curve truncated in order to rule out unrealistic demands . 

SUBROUTINE NORMAL* The normally distributed numbers are computed 
by the well known method of summing 12 U(0,1) random numbers, subtract^ 
ing the number six, and multiplying this quantity by the desired standard 
deviation and then adding the desired mean* The variable "FMEAN" and 
"STDDEV" are the mean and standard deviation desired for the distribution* 
Since it is possible to generate negative demands with this distribution 
the standard deviation and mean chosen should be such that the mean less 
twice the standard deviation is greater than zero* On the other hand, 
negative demands could be allowed to simulate materials returned to 
store or items to be exchanged. 

SUBROUTINE POISSON. The discrete analog of the integral method is 
used for this distribution. The probability mass function is summed 
term by term and compared to the U(0,1) random number obtained from the 
RAND subroutine. The Poisson number is taken to be the integer corresponds 
ing to that terra which when added to the sum of the previous terms ex= 
ceeds the U(0,1) random number. The variable "FLAMBDA" is the mean of 
the distribution. 
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SUBROUTINE NORMAL (IC0N6) 



C 

C 

COMMON INDX, NX S NRNC, NXIN, NRC, NSTRN 
C INSERT ALL THE DIMENSION AND COMMON CARDS OF THE MAIN PROGRAM 

GO TO (1,2), ICONS 
1 ITRAN = 1 
NRNC = lOOOO 
NRC = 0 
NSTRN = 0 
INDX = 6 
2 FMEAN = 50 o 0 
STDDEV = 25oO 

UPLIMIT = FMEAN + 2„*STDDEV 
FLWPT = FMEAN -2.*STDDEV 
DO 38 J = 1,12 
M = 13 +J 
Y = 0.0 
DO 5 L=l,12 
CALL RAND (ITRAN) 

RN = FLOATF( NXIN) /lOOOO. 

5 Y = Y + RN 

AO(M) = FMEAN + STDDEV*(Y“6. ) 

IF(AO(M) =FLWPT ) 30,30,35 
30 AO(M) = FLWPT 
GO TO 38 

35 IF(AO(M) - UPLIMIT) 38, 36, 36 
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36 AO(M) = UPLIMIT 
38 AO(M) = INTF( AO(M) + 0 5) 
RETURN 
END 



C 
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c 

c 

SUBROUTINE EXPONTL (ICON6) 

C 

C INSERT ALL THE DIMENSION AND COMMON CARDS OF THE MAIN PROGRAM 

COMMON INDX, NX, NRNC, NXIN, NRC, NSTRN 
GO TO (1,2), ICQN6 

1 ITRAN = 1 
NRC = 0 
NSTRN = 0 
INDX = 6 

2 NRNC = lOOOO 
THETA = .01 
ERROR = 0.01 

UPLIMIT = LOGF ( 1./ ERROR ) / THETA 
UPLIMIT = INTF( UPLIMIT ) 

DO 5 J=l,12 
CALL RAND (ITRAN) 

RN = FLOATF (NXIN) / 10000. 

M = 13 + J 

AO(M) = INTF( L0GF(1./(1. - RN))/THETA +0.5) 

5 AO(M) = MIN1F ( AO(M) , UPLIMIT ) 

RETURN 

END 
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c 

c 

SUBROUTINE POISSON (ICON 6) 

C 

C INSERT ALL THE DIMENSION AND COMMON CARDS OF THE MAIN PROGRAM 

COMMON INDX,NX,NRNC,NXIN,NRC,NSTRN 
GO TO (1,2), IC0N6 

1 ITRAN = 1 
INDX = 6 
NSTRN = 0 
NRC = 0 

2 NRNC = lOOOO 
FLAMBDA = 4. 

IUPLIM = 12 

DO 11 J = 1,12 
CALL RAND (ITRAN) 

RN = FLOAT F (NXIN) / 10000 „ 

TESTNR =RN tt EXPF( FLAMBDA) 

M=13+J 

IF(TESTNR-1.)5,5,8 
5 A0(M)=0.0 
GO TO 11 
8 PRODUCTS. 0 
Y=1.0 

DO 9 K = 1, IUPLIM 
XNR1 = FLOATF (K) 

FRODUCT = PRODUCT »XNR1 
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Y = FLAMBDA“*K/FR0DUCT + Y 
IF ( TESTKR-Y ) 10 , 10 , 9 
9 CONTINUE 

10 AO(M) =FLOATF( K ) 

11 CONTINUE 
RETURN 
END 
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c 

c 

SUBROUTINE UNIFORM (ICON6) 

C 

C COMMON STORAGE 

C INSERT ALL THE DIMENSION AND COMMON CARDS OF THE MAIN PROGRAM 

COMMON INDX, NX, NRNC, NXIN, NRC , NSTRN 
GO TO (1,2), ICON6 

1 ITRAN = 1 
NRC = 0 
NSTRN = 0 
INDX = 6 

2 NRNC = lOOOO 
FLWFT =50.0 
UPENDPT = 80.0 
DO 5 J = 1,12 
CALL RAND (ITRAN) 

RN = FLOATF (NXIN) / 10000, 

M = 13 + J 

5 AO(M) = INTF( FLWPT +( UPEND FT - FLWPT)*RN + 0.5) 

RETURN 

END 

C 
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c 

c 

SUBROUTINE RAND (ITRAN) 

C 

C RANDOM NUMBER GENERATOR 

C 

C LOCAL STORAGE 

DIMENSION XI (8) ,NXI(8) 

C INSERT ALL THE DIMENSION AND COMMON CARDS OF THE MAIN PROGRAM 

COMMON INDX,NX,NRNC,NXIN,NRC,NSTRN 
EQUIVALENCE (XI ,NXI) , (T1,NT1) , (INDX,DX) , (D1,ND1) 

C 

GO TO (10, 20), ITRAN 
10 ITRAN=2 

B Xl(l)=3514 6016 2524 6131 

B XI(2)=0337 1363 2712 7740 

B XI(3)=1760 3011 0710 3016 

B XI(4)=0670 1154 5656 1316 

B XI(5)=1506 7663 7414 1566 

B Xl(6)=2100 7160 3140 2631 

B XI(7)=1037 4327 7340 4007 

C Xl(8)=0 

C 

20 NRC=NRC+1 
L=INDX+1 
ND1=L 

B DX=D1*7 
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M=INDX+1 



IID1=M 

B T1=D1*7 

N=NT1+1 

NXI(M)=NXI(L)+NXI(N) 

B XI(M)=XI(M)*3777 7777 7777 7777 

NX=NXI(M) 

30 lda(nx) ,muf(nrnc) ,sta(nxin) . 

c 

RETURN 

END 

C 



105 



TAELE 1 



RESULTS 


FROM TEST RUNS 


ON THE RANDOM 


DEMAND SUBROUTINES 




NORMAL DISTRIBUTION 




SAMPLE SIZE 500 




MEAN 


XBAR 


STDDEV 


SAMPLE VAR 


Vs. V. 


25 


26,54 


25 


550.1 


23.4 


50 


49 *44 


25 


615.6 


24.8 


75 


74.79 


25 


561.3 


23.7 


100 


101,71 


25 


547.7 


23.4 


125 


124.34 


25 


506.99 


22,5 


150 


149.50 


25 


511.69 


22.6 


EXPONENTIAL DISTRIBUTION 




SAMPLE SIZE 500 




MEAN 


XBAR 


STDDEV 


SAMPLE VAR 


Vs. V. 


100 


102.6 


100 


9,657.9 


98.30 


20 


19.78 


20 


357.82 


18.93 


10 


9.89 


10 


88.86 


9.42 


4 


4.04 


4 


15.16 


3.90 


2 


1.88 


2 


3.46 


1.83 


POISSON DISTRIBUTION 




SAMPLE SIZE 500 




MEAN 


XBAR 


STDDEV 


SAMPLE VAR 


Vs. V. 


4 


4.10 


2.00 


3.99 


2.00 


5 


5.06 


2.24 


5.27 


2.29 


6 


5.92 


2.45 


5.10 


2.26 


7 


7.10 


2.64 


7.15 


2.67 


8 


7.83 


2.83 


7.51 


2.74 


9 


9.06 


3.00 


9.54 


3.09 


10 


9.98 


3.16 


9.74 


3.12 
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UNIFORM DISTRIBUTION 



INTERVAL 


IE AH 


XBAR 


(0,100) 


50 


50.91 


(10,90) 


50 


49.81 


(20,80) 


50 


49.07 


(30,70) 


50 


49.74 


(40,60) 


50 


49.13 



SAMPLE SIZE 500 



STDDEV 


SAMPLE VAR 


C 

< 


28.86 


823.59 


28.69 


23.08 


526.92 


22,95 


17.32 


273.34 


16.53 


11.53 


135.28 


11.62 


5.77 


32.13 


5.66 
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APPENDIX II 



DICTIONARY OF VARIABLES 



Variable 






Name 




Definition 


Unit 


A 




Production setup cost 


Dollars/Setup 


AO 




Input data array 


None 


APB 




Branch physical storage cost 


Cents/Product unit/ 








Year 


APP 




Whs physical storage cost 


Cents/Product unit/ 








Year 


AZR 




Opportunity cost 


Cents/Dollar/Year 


BASE 




Product code number 


None 


BCCJK 


J,K 


Broken case count, item-location 


Unit 


BCK 


K 


Branch code 


None 


BCP 




Broken ctn and pallets, system 


Unit 


BCFCJ 


J 


Broken ctn and pallet count, item 


Unit 


BCPCJK 


J,K 


Broken ctn and pallet count 


Unit 


BCPI 




Broken ctn and pallets, product 


Product Unit 


BPCJK 


J,K 


Broken pallet count, item-location 


Unit 


bpijk 




Computational factor 


None 


BSI 




Batch size increment 


Product Unit 


BSM 




Batch size, minimum 


Product Unit 


BF 


J,K 


Branch usage flag 


None 


BG 




Branch requisitioning goal 


Days (of Supply) 


BPCI 




Broken pallet count, product 


Unit 


CBC 




Cost of broken ctns, system 


Dollars 


CBCC 




Cost of picking broken ctns 


Dollars/Unit 


CBC I 




Cost of broken ctns, product 


Dollars 


CBCJ 




Cost of broken ctns, item 


Dollars 


CBCJK 




Cost of broken ctns, item-location 


Dollars 


CBP 




Cost of broken pallets, system 


Dollars 


CBPC 




Cost of picking broken pallets 


Dollars/Unit 


CBPI 




Cost of broken pallets, product 


Dollars 


CBPJ 




Cost of broken pallets, item 


Dollars 


CBPJK 




Cost of broken pallets, item-loca- 








tion 


Dollars 


CBO 




Cost of branch orders, system 


Dollars 


CBOK 




Cost of branch orders, branch 


Dollars 


CF 




Freight cost, system 


Dollars 


CFK 


K 


Freight cost, branch 


Dollars/100 Pounds 


CH 




Holding costs, system 


Dollars/Year 


CHI 




Holding costs, product 


Dollars/Year 


CHJ 




Holding costs, item 


Dollars/Year 


CHJJ 




Opportunity plus physical holding 


Dollars/Unit/Week 


CHJK 




Holding costs, item-location 


Dollars/Unit/Week 


CHK 


K 


Holding cost rate, branch 


Dollar s/Unit/Week 


CIJ1 




Whse unit cost 


Dollars/Unit 


CIJ 


K 


Branch unit cost 


Dollars/Unit 


CJ 


J 


Production unit cost 


Dollars/Unit 


CL 




Cost of losses, system 


Dollars 


CLI 




Cost of losses, product 


Dollars 
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CU 




Cost of losses, item 


Dollars 


CUK 




Cost of losses, item-location 


Dollars 


CLK 


K 


Cost of losses, branch total 


Dollars 


CKZJK 




Variance factor 


None 


CM2 




Order rate multiplier 


None 


CP 




Cost of setups, system 


Dollars 


CPD 




Cost of loss per demand, system 


Dollars/Unit 


CPDK 




Cost of loss per demand, branch 


Dollars/Unit 


CPP 




Cost of order preparation and 








process, system 


Dollars 


CPPK 




Cost of order preparation and 








process, branch 


Dollars 


CPK 


K 


Processing cost, branch 


Dollars 


CPFP 




Cost of picking full pallet, rate 


Dollars/Pallet 


CPLK 


K 


Cost of picking and loading, branch 


Dollars 


CPPJ 




Cartons per pallet 


Ctns/Pallet 


CPFO 




Order preparation and processing 


Dollars/Order 






cost, rate 


CPIJK 


J,K 


Processing cost, item-location 


Dollars 


CR 




Cost of receiving (processing) per 


Dollars/lOO Pounds 






unit 


CSIJK 




Profit (selling-branch cost) also 


Dollars/Unit 






cost of shortages for both direct 
shortages at branches and short- 
ages resulting from stockout at 
the central warehouse 




CTOT 




Sum of CPL, CP, CL, SH, SYS-TOT 


Dollars 


CSZ 


J 


Cost of losses, item 


Dollars 


CPL 




Cost of picking, loading, system 


Dollars 


DBI 




Demand rate, branches less warehouse 


Product unit/Week 


UBIJJ 




Demand rate, branches less warehouse 


Unit /Week 


Dt 




Predicted demand, product 


Product unit/Week 


DIJ 


J 


Predicted demand, item 


Unit/Week 


DIJK 


J,K 


Predicted demand, item-location 


Unit/Week 


DICHI 




Mean holding cost/unit of product 


Dollars/Product unit 






times system demands (DI) 




DIFR 




Mean demand between orders, product 


Unit 


DIFRFR 




Mean demand between orders, item 


Unit 


DLBI 




Leadtime demands, branches less 


Product unit 






whse, product 




DLBJ 




Leadtime demands, branches less 


Unit 






whse, item 




DLI 




Leadtime demands, whse, product 


Product unit 


DLIJ 




Leadtime demands, whse, item 


Unit 


DM 


J,K,L 


Monthly demands 


Unit 


DIR 


J 


Probable demand in half of leadtime 


Unit 


DSCAJ 




Discrete adjustment interval 


Week 


DSIJ 




Satisfied demand, item 


Unit 


DSIJK 


J,K 


Satisfied demand, item-location 


Unit 


DWIJK 


J,K 


Weekly demands, item-location 


Unit 


DZ 




Sum of DZIJK over all products 


Unit 


DZIJK 




Mean of distribution for number of 


None 



demands in prodcn cycle plus 
leadtime 
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F 




Correction factor to Wilson's EOQ; 
Sum of all FIJK 


Product unit 


F2 




Square of F 


Product unit squared 


FITEM 


J 


Item number, floating point mode 


None 


FIJK 




The demands in a production lead- 
time adjusted for the variances 
of the demand and leadtime dis- 
tributions 


Unit 


FKL 




Multiplier, branch requisition 
objective floor 


None 


FK2 




Multiplier, branch reorder level 
floor 


None 


FK3 




Multiplier, warehouse requisition 
objective floor 


None 


FK4 




Multiplier, warehouse reorder level 
floor 


None 


FK5 




Multiplier, item production trigger 
level 


None 


FK 6 




Multiplier, production trigger 
level, warehouse 


None 


FKIJK 




Variance factor for demands and 
production reorder leadtime j 
used to compute RIJK(J,K) for 
branches other than warehouse 


None 


FKZJK 




Same as FKIJK except applies to 
branch interim reorder leadtime j 
used to compute RIJK(J,l) 


None 


FL 




Warehouse production leadtime 5 for 
central whse LC + LP + IRl/2 


Week 


FLC 




Communication leadtime 


Week 


FLIK 


K 


Branch production leadtime ; Sum of 
LRl/2,LC,LP,LSK,LW/2 


Week 


FLDI 




Maximum EOQ due to shelf life 
limitations 


Product unit 


FLP 




Production-run leadtime 


Week 


FLPDTM 




Time of placing most recently de- 
livered production allocation 


A particular week 
Unit ^ 


FLQI 




Optimal ratio of EOQ to square root 
of product demands 


FIR 




Reaction leadtime 


Week 


FLRI 




Review interval 


Week 


FLKIJK 




Variance correction factor for 
demands and production leadtime 
to be applied to production lead- 
time demands 


None 


FLZK 


K 


Branch interim reorder leadtime 5 
LSK + LRI/2 + LW/2 


Week 


FLS 




Leadtime, shipping, mean 


Week 


FLSK 


K 


In-transit leadtime, branch 


Days 


FMU 




System capacity constraint 


None 


FMUK 


K 


Branch capacity constraint 


None 


FMDF 




Demand forecast multiplier 


None 


FMBG 




Branch goals, multiplier 


None 


FMPG 




Plant goals, multiplier 


None 
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FMCJ 




Standard cost multiplier 


None 


FMLP 




Leadtime multiplier, plant to 
warehouse 


None 


FMLWB 




Leadtime multiplier, warehouse to 
branch 


None 


FNOPR 




Number of production runs, floating 
point 

Freight rate 


None 


FRK 


K 


Cents/Cwt 


FZJK 




Demands in branch interim order 
leadtime 


Unit 


GENDEM 




Subroutine to generate demands 


None 


HC 




Holding cost due to rounding to 
pallets 


Dollars 


ICON 




Parameter for RXCOKP subroutine 


None 


icon 




Rule selector; 1-SRI, 2-other 


None 


ICON2 




Not used 


None 


ICOK3 




Stockout policy selector, 

1- issue to zero balance 

2- do not make partial issues 


None 


IC0N4 




Round to cartons selector, 

1- do not round to full case lot 

2- round to next higher case lot 


None 


ICON 5 




Round to pallets selector, 

1- do not 

2- round if holding cost is less 
than handling costs 


None 


I 




Index 


None 


IE OF 




End of file, fixed point mode 


None 


IT 




Index 


None 


J 




Index for item 


None 


J1 




Index used in subroutine predict 


None 


J2 




Index used in subroutine predict 


None 


K 




Index for branch 


None 


L 




Index for month and production order 


None 


M 




Index and parameter for subroutine 
QUEUE P 


None 


N 




Index parameter for subroutine QUEUEB None 


NJ 




Number of items 


None 


NK 




Number of branches 


None 


OC 




Order count, system 


Orders 


XI 




Order count, product 


Orders 


XJ 




Order count, item 


Orders 


XJK 


J,K 


Order count, item-location 


Orders 


ODD 




Order delivery date 


Week 


OHJK 


J ,K 


On-hand, item-location 


Unit 


FSK 


K 


Shipping frequency 


Week 


OOJK 


J,K 


On order, item-location, cumulative 


Unit 


ORD 




Order quantity 


Unit 


XK 




Order count, branch, mean 


Orders 


PAGE 




Page number 


None 


PBO 




Subroutine to place branch orders 


None 


FG 




Plant goal 


Days (of Supply) 


PJ 


J 


Unit selling price, branch 


Dollars/Unit 


P0A3 


J 


Production order allocation to item 


Unit 
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POAW 


J 


Production order allocation to 
warehouse 


Unit 


FOBS 


J 


Production batch size 


Product unit 


?ODD 




Production order delivery date 


Week 


rc 




Picking cost due to roundg to 
cartons 


Dollars 


PRIJK 




Probability of order between produc- 
tion allocations 


None 


PREDICT 




Subroutine to predict demands 


None 


QI 




Economic prodcn order qty 


Product unit 


QIZ 




Computational factor for QI 


Product unit 


QIZZ 




Minimum of alternate values of EOQ 


Product unit 


QIFR 




Computational factor for RIJK, XIJK 


Product unit 


QlfR 




The number of orders corresponding 
to the adjusted EOQ quantity 


None 


QJ 


I 


Order item index in subroutine queue 


None 


QK 


I 


Order branch index in subroutine 
queue 


None 


QQRD 


I 


Order quantity in subroutine queue 


Unit 


QODD 


I 


Order delivery date in subroutine 
queue 


A specific week 


QUEUEB 


N 


Subroutine for queuing branch orders 


None 


QPOAS 


J,I 


Production item allocatn in queue 


Unit 


Q POAW 


J,I 


Production whse alloctn in queue 


Unit 


QPOBS 


I 


Production batch size in queue 


Product unit 


QPODD 


I 


Production order date in queue 


A specific week 


QUEUE P 


M 


Subroutine for queuing production 
order 


None 


QRIJK 




Modified EOQ value 


Product unit 


R 




Scaling factor j POAS/CRD 


None 


RBI 




Requisition objective level, branch 
only 


Product unit 


RGX 




Ratio for alternate trigger level 


None 


RI 




Requisition objective level (goal), 
product 


Product unit 


RIJ 




Requisition objective level (goal), 
item 


Unit 


RIJJ 




Requisition objective level (goal), 
item 


Unit 


RIJK 


J,K 


Requisition objective level (goal), 
item-location 


Unit 


RJK 




Maximum of RZZJK and XIJK plus lead- 
time demands 


Unit 


RJ1 




Reauisition objective level, ware- 
house only, item 


,Unit 


RLT 




Ratio of losses to total demands, 
sys tem 


None 


RLTI 




Ratio of losses to total demands, 
product 


None 


RLTJ 




Ratio of losses to total demands, 
item 


None 


RLTJK 




Ratio of losses to total demands, 
item-location 


None 
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RPJK 




Computational factor for requisition 


Unit 








objective 






RST 




Stockturn ratio, system 


None 




RSTI 




Stockturn ratio, product 


None 




RSTJ 




Stockturn ratio, item 


None 




RSTK 




Stockturn ratio, branch 


None 




RSTJK 




Stockturn ratio, item-location 


None 




RUNNO 




Run number 


None 




RXCOMP 




Subroutine to compute SRI rules 


None 




RZJK 




Alternate value of requisition 


Unit 








objective 






RZZJK 




Minimum of alternate values of RIJK 


Unit 




RDTP 




Subroutine to read input data 


None 




SC 


J 


Size (item) code 


None 




SD 




System demands 


Product 


unit 


SDB 




System demds less warehouse demds 


Product 


unit 


SDL 




System demds in leadtime 


Product 


unit 


SDLB 




Sys demds in leadtime, branches only 


Product 


unit 


SG 




System goal all products, and 


Product 


unit 






locations 






SGP 




System goal, less warehouse 


Product 


unit 


SGPP 




System goal, warehouse 


Product 


unit 


SI 




Inventory position, product 


Product 


unit 


SIJ 


J 


Inventory position, item 


Unit 




SIJK 


J,K 


Inventory position, item-location 


Unit 




SJ 


J 


Product unit conversion factor 


Product 


unit/Unit 


soc 




Stockout count, system 


None 




SOCI 




Stockout count, product 


None 




SOCJ 




Stockout count, item 


None 




SOCJK 


J,K 


Stockout count, item-location 


None 




SGRD 




Cumulative ord qty 


Product 


unit 


SOW 




Stockout weeks, system 


Week 




SOWI 




Stockout weeks, product 


Week 




SOWJ 




Stockout weeks, item 


Week 




SOWJK 


J,K 


Stockout weeks, item-location 


Week 




SQ 




Economic order quantity, system 


Product 


unit 


SQI 




Economic order quantity, product 


Product 


unit 


SQP 




Economic order quantity, branches 


Product 


unit 






only 






SQWI2 




Square of Wilson's EOQ, product 


Product 


unit 


SR 




System requisition objective, all 


Product 


unit 






prods + Iocs 






SRP 




System requisition objective, 


Product 


unit 






branches only 






STI 




Coverage factor, product 


Week 




STIK 




Coverage factor, location 


Week 




SY 




Average inventory position, system 


Product 


unit 


SYI 




Average inventory position, product 


Product 


unit 


SYJ 




Average inventory position, item 


Unit 




SYJK 


J,K 


Average inventory position, item- 


Unit 








location 






SOWK 


K 


Stockout weeks, branch 


Week 




T 




Time index 


None 
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TI 




Coverage provided by QI (BAR) plus 


Week 






one half of reaction time 




TP 




Temporary storage location for OHJK 


None 


TP2 




Temporary storage location for DSIJK 


None 


TX 




Shelf life 


Month 


T12 




Interval between time of placing 


Week 






most recently delivered production 








allocation and present time 




UL 




Units loss, system 


Product unit 


ULI 




Units loss, product 


Product unit 


UU 




Units loss, item 


Unit 


ULJK 




Units loss, item-location 


Units 


UPCJ 


J 


Units per carton 


Units/Carton 


VARDF 




Function to compute variance 


None 


VDIJ1 




Rel-variance of demands, whs 


None 


VDFR 


J 


Probability of interim order 


None 


VPIJK 




Rel-variance of demand, item-location 


None 


VDS 




Value of demand satisfied, system 


Dollars 


VDSI 




Value of demand satisfied, product 


Dollars 


VDSJ 




Value of demand satisfied, item 


Dollars 


VDSK 


K 


Value of demand satisfied, branch 


Dollars 


VDSJK 




Value of demand satisfied, item- 


Dollars 






location 




VLP 




Rel-variance of prodcn leadtime 


None 


VLZK 




Rel-variance of branch leadtime 


None 


VSIT 




Value of inventory in-transit 


Dollars 


VSYI 




Value of ave inventory, product 


Dollars 


VSYJ 




Value of ave inventory, item 


Dollars 


VSYK 


K 


Value of ave inventory, branch 


Dollars 


VSYJK 




Value of ave inventory, item-location 


Dollars 


VSY 




Value of ave inventory, system 


Dollars 


VY3IT 




Value of ave inventory in-transit 


Dollars 


VZD 




Value of total demands, system 


Dollars 


VZDI 




Value of total demands, product 


Dollars 


VZDJ 




Value of total demands, item 


Dollars 


VZDK 


K 


Value of total demands, branch 


Dollars 


VLIK 


K 


Rel-Variance of branch leadtime 


None 


VDZJK 




Rel-variance of demands, item- 


None 






location 




V2 




Control knob - system requisition 


None 






level in weeks of demand except 








for central whse under SRI rules 




V3 




Control knob - system requisition 


None 






level in weeks of demand includ- 








ing central wh 3 e under SRI rules 




V4 




Control knob - difference between 


None 






V2 and W1 




V 5 




Control knob - difference between 


None 






V3 and W2 




V8 




Control knob - expected maximum on 


None 






hand level, all branches 




V9 




Control knob - expected average on 


None 



hand level, all branches 
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V10 




Control knob - same as V8 except 
warehouse is excluded 


None 


Vll 




Control knob - same as V9 except 
warehouse is excluded 


None 


w 




Temp storage location 


None 


WJ 


J 


Unit weight 


Pounds /Unit 


WPSO 




Weeks per stockout, system 


Week 


WPSOI 




Weeks per stockout, product 


Week 


WPSOJ 




Weeks per stockout, item 


Week 


WPSOJK 




Weeks per stockout, item-location 


Week 


W1 




System requisition level in weeks 
of demand except for central 
whse under alternate rules 


Week 


W2 




System requisition level in weeks 
of demand including central 
whse under alternate rules 


Week 


X 




Temporary storage location 


None 


XCI 


J 


Dealer const (unit plus shipping 
plu3 receiving) 


Dollars/Unit 


XDIJ 


J 


Demands over shelf life, item 


Unit 


XDIJK 


J,K 


Demands over shelf life, item-loca- 
tion 


Unit 


XI 




Reorder trigger level, product 


Product unit 


XIJ 


J 


Reorder trigger level, item 


Unit 


XIJK 


J,K 


Reorder trigger level, item-location 


Unit 


XJ1 


J 


Reorder trigger level, warehouse 


Unit 


XJK 




Reorder trigger level, item-location 
other than warehouse 


Unit 


XK2IK 


J 


Interim order leadtime times reorder 
floor multiplier (FK2) 


None 


XRIJK 


J,K 


Requisition objective level at a 
branch in weeks of item location 
demands 


Week 


XRIJZ 


J 


Requisition objective level at ware- 
house in weeks of total item de- 
mand 


Week 


XXIJ 


J 


Reorder level for item in weeks of 
total item demand 


Week 


XXIJZ 


J 


Reorder level for warehouse in weeks 
of total item demand 


Week 


Y 




Temporary storage location 


None 


ZDIJ 




Total demands item 


Unit 


ZDIJK 


J,K 


Total demands item-location 


Unit 


ZVC 




Total of CH, 0C, CL; system 


Dollars 


ZVCK 




Total of CHK, CBOK, CLK; branch 


Dollars 


Z 




Obsolescence cost (not used in 
present form of simulator) 


Cents/Unit 



140 









*- 





Approved for public release y 
distribution unlimited.. 



